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Abstract 

We study the solutions of linear Schrodinger equations in which the potential energy 
is a periodic function of time and is sufficiently localized in space. We consider the po- 
tential to be close to one that is time periodic and yet explicitly solvable. A large family 
of such potentials has been constructed and the corresponding Schrodinger equation 



solved by Miller and Akhmediev [22|. Exact bound states, or breather modes, exist in 
the unperturbed problem and are found to be generically metastable in the presence 
of small periodic perturbations. Thus, these states are long-lived but eventually decay. 
On a time scale of order e -2 , where e is a measure of the perturbation size, the decay is 
exponential, with a rate of decay given by an analogue of Fermi's golden rule. For times 
of order e _1 the breather modes are frequency shifted. This behavior is derived first 
by classical multiple scale expansions, and then in certain circumstances we are able 



to apply the rigorous theory developed by Soffer and Weinstein [28] and extended by 
Kirr and Weinstein |15| to justify the expansions and also provide longer time asymp- 
totics that indicate eventual dispersive decay of the bound states with behavior that 
is algebraic in time. As an application, we use our techniques to study the frequency 
dependence of the guidance properties of certain optical waveguides. We supplement 
our results with numerical experiments. 
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1 Introduction and Overview 



We are interested in the initial value problem for the linear Schrodinger equation in one 
space dimension 



U),f= ( -1% + V{x,t) )/ = Hit) f (1.1) 



Here, / is a complex- valued function ofx 6 R and t e R. We assume that V(x, t) is a smooth 
real-valued potential energy function which is sufficiently localized in space (for example, of 
Schwartz class). In our specific applications, V(x,t) will be taken to be a periodic function 
of t. However, the techniques we use can be adapted for to more general time-dependence 



15 1 . Note that ( |1 . 1| ) is a nonautonomous Hamiltonian system: 



id tf = (1-2) 



where h denotes the Hamiltonian energy: 
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h[f,f,t} = j l^\d x f\ 2 + V(x,t)\f\ 2 j dx. (1.3) 

If V(x, t) is not independent of t, h is not a conserved integral of the flow. On the other hand, 
since the potential V is real- valued, the flow defined by ( |1.1|) is always unitary in L 2 (IR), i.e. 

J \f{x,t)\ 2 dx = J |/(x,0)| 2 dx, teM. (1.4) 

In applications, it is often natural to decompose V(x, t) as: 

V(x,t) = V (x,t) + W(x,t), 

where V (x,t) denotes an unperturbed potential, and W(x,t) = V(x,t) — V (x,t), denotes a 
small perturbation. Thus, 

H{t) = H {t) + W(t), (1.5) 



and ( |1.1| ) can be rewritten as: 

idtf = (H (t)+W{t) )/ (1.6) 

Here, we have denoted the multiplication operator / > W(x,t)f by W(t). The choice of 
V (x,t) is often dictated by some a priori knowledge of the solutions of the unperturbed 
system 

idtf = H (t)f. (1.7) 



A problem of importance is then to contrast the detailed dynamics of solutions to (|1.6| ) 



with those of the unperturbed system ( |1.7| ). In particular, if (\l.lQ has bound state solutions 
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(breather modes, or solutions having finite energy and not decaying as \t\ — > oo) do they 
persist in the perturbed dynamical system (\1.6J ? 

The simplest variant of this problem is the case where the unperturbed part is stationary, 
i.e. Vo(x, t) = Vo(x). Suppose the operator TCo has an L 2 eigenf unction. The unitary 
evolution of the spatially localized eigenfunction is time-periodic and represents a bound 



state solution of the unperturbed Schrodinger equation ( |1.7|) . The perturbed model (in this 



and in the more general case when Ti, has multiple discrete eigenvalues) is related to the 
problem of ionization of an atom by a time-dependent electromagnetic field |L7|, [Hj and the 
problem of describing the effects of weak inhomogeneities on the propagation of continuous 



waves in optical fibers p0 |. Using a time-dependent method developed in the context of (i) 



quantum resonances and the perturbation theory of embedded eigenvalues in the continuous 



spectrum [26] and (ii) resonances and radiation damping of bound nonlinear wave equations 



27j, Soffer and Weinstein studied the metastability of such states p8|. Generalizations of 



this theory for handling multifrequency perturbations |T5| and the interference of multiple 
bound states in the unperturbed problem Jl6| have been explored by Kirr and Weinstein. 
Based on the observation that the mechanism for instability of the bound state is coupling 
of the bound state to the continuous spectral modes, the analysis was carried out at the 
level of the coupled equations for the bound state and dispersive components of the solution. 
Under general hypotheses on unperturbed Hamiltonian (local energy decay estimates on the 
unitary propagator e~ lt °) this equivalent dynamical system was studied and it was shown 
that a bound state is generically unstable but long-lived. The lifetime is given by a formula 
analogous to the Fermi golden rule 



In this paper we consider the case where the unperturbed Hamiltonian is genuinely time- 
dependent. A physical application of the theory we develop, in the context of frequency 
detuning in periodically modulated optical waveguides 0, [|, will be presented in §|5|. Let 
Ho(t) = —\d"l + V (x, t), where V (x, t) is smooth, periodic in t with the same period for 
each x and of sufficiently rapid decay for large x for each i £ 1. The particular choices of 
Vo(x, t) we consider in this paper belong to a large family of very special, so-called separable, 
time-dependent potentials, Vo(x,t), x G R, studied by Miller and Akhmediev |22| . The 
separable potential Vq(x, t) can be chosen to be time-periodic, in which case the unperturbed 
problem supports exact bound states (breather modes) and the initial value problem for flLTj ) 
can be solved exactly. That is, a complete set of eigenmodes and generalized eigenmodes 
can be explicitly displayed with respect to which the dynamics of ( |1.7|) is diagonal. This 
class of potentials is intimately connected with the soliton theory of completely integrable 
multicomponent cubic nonlinear Schrodinger equations [19|. |9|] . 

The existence of such exact breather modes in the unperturbed time-periodic problem 
is quite remarkable and we believe that this is a highly non-generic phenomenon^. Indeed, 
from a general dynamical systems perspective, (|1.7|) with such a choice of Vo(x,t), may be 
viewed as a parametrically forced wave equation (here we are actually considering the time- 



1 The scarcity of breather solutions of nonlinear wave equations defined on a spatial continuum of infinite 
extent has been extensively explored in the setting of perturbations of the completely integrable sine-Gordon 
equation; see, for example, || [|, ||, H. The connection with linear nonautonomous problems can be made 
by viewing a breather solution of a nonlinear dynamical problem as a bound state of a linear problem with 
a given (self-consistent) potential. 
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periodic function V (x, t) itself to be the sum of a time-independent part and a time-periodic 
modulation). One therefore expects that the presence of resonances will perturb the Floquet 
multipliers (corresponding to bound states) off the unit circle as in the elementary example 
of Mathieu's equation [0. The persistence of breather solutions under the time-periodic 
perturbation would imply the non-departure of a Floquet multiplier from the unit circle to 
all orders in the size of the perturbation. The fact that infinitely many such conditions hold 
for these special separable potentials is no doubt linked to the infinite sequence of symmetries 
and time invariants enjoyed by the completely integrable nonlinear flow that underpins the 
construction of the separable potentials (see |Appendix A| for more details). Of course, this 
is only a heuristic picture. In fact, the perturbation theory of the Floquet multipliers is 
complicated by the fact that they are embedded in the continuous spectrum which covers 
the unit circle. However, spectral deformation methods have been developed for some classes 
of models that could well be adapted here. Relevant technical details can be found in || |TTl|p| 

We want to make our motivation for pursuing deformations of these admittedly rather 
special periodic potentials very clear. First of all, the problem is relevant to the analy- 
sis of optical waveguides. In the paraxial approximation, the slowly varying envelope of a 
highly oscillatory electric field in a dielectric medium with inhomogeneous dielectric proper- 
ties (index of refraction) satisfies an equation of form ( [Lip . Here, t denotes the longitudinal 
variable, the direction of propagation, and x is the transverse spatial variable. For inho- 
mogeneous index profiles corresponding to exactly solvable potentials, light of a particular 
frequency propagates as a non-attenuating bound state mode in these wave guides. However, 
if the light frequency deviates from the "integrable frequency" the propagating wave will be 
governed by the perturbed equation (|1.6| ). Thus the question of whether such modes persist 
and if not what their lifetime is for the perturbed dynamics naturally arises. We will give 
more details about this problem in §[5|. 

But it is also true that the study of perturbed separable periodic potentials is important 
in general terms. Given an arbitrary time-periodic potential in the Schrodinger equation, 
one wants to study the corresponding dynamics using perturbation theory. In doing so, the 
first question one must address is that of finding a "nearby" problem that can be solved 
exactly. We simply take the point of view that many periodic potentials will be closer 
to a separable periodic potential (in a sense that can be made precise) than to any time- 
independent potential. 

In any case, with the explicit spectral theory associated with V (x, t) in hand, our goal 
is to carry out a detailed analytical study of the coupled mode dynamics induced by a time- 
dependent perturbation W(x,t). We establish the generic metastabilty of the exact bound 
states associated with separable periodic potentials Vq(x, t) and obtain a detailed picture of 
the dynamics. 

The paper is structured as follows. In we first review the construction of time- 

2 For nonlinear wave equations defined on an infinite lattice (e.g. discrete sine-Gordon, discrete </> 4 ), 
breather solutions can be constructed for sufficiently large lattice spacing; see for example . The radiative 
decay of such discrete breathers, for sufficiently small lattice spacing, is expected to be governed by a 
mechanism of the kind studied in this paper; see also |2?|] . Related to this are results for the dynamics of 
kinks of discrete nonlinear wave equations, in which the techniques of this paper have been used to study the 
"pinning" of discrete kinks on a lattice site. This pinning is marked by the slow radiative decay of spatially 
localized and time-periodic or quasiperiodic oscillations about a static kink [Of. 
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dependent exactly solvable potentials from a set of discrete data, and then show how 



the initial value problem for such separable potentials can be solved explicitly. We then 
describe how properties of the separable potentials depend on the choice of the discrete data 
generating them. Next, we derive by projection onto an orthonormal basis the general cou- 
pled mode equations which arise when a separable potential is perturbed by some arbitrary 
correction W(x,t). This section will then conclude with a detailed derivation of the two- 
soliton time-periodic potential and its associated explicit complete set of bound states and 
generalized eigenfunctions. More details about the separable potentials described in §|2] are 
given in |Appendix Aj . 

In §Q Floquet theory is then used to map the coupled mode equations to a system asso- 
ciated with a time-dependent perturbation of an autonomous system, a situation analyzed 
in detail in |28| and In §[|, we then describe the dynamics of solutions of the coupled 



mode equations for even time-periodic perturbations W(x,t) of a separable two-soliton even 
time-periodic potential Vo(x,t) (Vo(x,t) and W(x,t) both share even parity in x and have 
the same temporal period). In particular, we study the initial value problem when the ini- 
tial condition is a pure bound state of the unperturbed problem. First, we study the small 
time behavior of the coupled mode equations (without requiring W(x,t) to be small) and 
deduce that the bound state amplitude behaves as 1 — Ct 2 for some constant C and interpret 
this result in the context of the theory of ideal measurements in quantum mechanics (the 
"watched pot" effect). We then assume the perturbation W(x,t) to be small, of size e, and 
seek the behavior of the bound state amplitude over intermediate times of order e -1 and e~ 2 
using the classical method of multiple scales. 

We show the existence of a perturbation-induced frequency shift of the breather mode 
evident on timescales of order e -1 and exponential decay of the bound state mode amplitude 
on timescales of order e~ 2 . The condition for the decay constant to be nonzero is a direct 
analog of the "Fermi golden rule" . 

Then, using the transformation to an autonomous system found in §|3|, we show how the 
rigorous theory developed for multifrequency perturbations of autonomous systems by Kirr 
and Weinstein |15j can be applied in some cases to justify the multiple scales calculation, 
and to provide more detailed information about the infinite time behavior of the solution. 
This analysis completes the portrait of the dynamics, showing that the exponential decay 
is ultimately washed out in a sea of dispersive waves, at which point the decay becomes 
algebraic in time. 

Having described the theory, in §[| we consider an application of the analysis to a problem 
of frequency detuning in planar optical waveguides. Finally, in §|6] the prediction of an 
exponential decay constant T for the bound state mode amplitude found in section f| is 
compared to numerical simulations of the perturbed time- dependent Schrodinger equation. 

A detailed description of the theory of separable potentials, at once summarizing for 
completeness and also further developing the results of [53], can be found in [Appendix A. 



In |Appendix B| the reader will find the proofs of the decay estimates that we will use in 

in order to apply the results of [Hj]. 

Regarding notation: We will use the inner product 

POO 

(/(■),»(•)>= / f{xy g (x)dx (i. 
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on L 2 (R). Occasionally, the angled brackets will denote the inner product in more gen- 
eral Hilbert spaces. Linear operators will be denoted with calligraphic letters, vectors with 
arrows, and matrices with boldface letters. We will often use the function defined by: 

(x) = (l+x 2 ) 1/2 . (1.9) 

Complex conjugation will be denoted with stars, and time averages will be denoted with 
bars. 



2 Exactly solvable time-dependent potentials 

In this section we recall for our purposes a class of time-dependent potentials V (x, t) related 
to M-soliton solutions of certain completely integrable nonlinear flows. Because of the inti- 
mate connection of these potentials to integrable systems, it is possible to explicitly derive 
the spectral representation associated with such potentials |22j . This section is divided into 



five parts. First, the direct construction of separable potentials from a set of discrete data 
T> is outlined. Then, we show how the same discrete data T> gives rise to formulas for a 
complete set of modes for the time-dependent Schrodinger equation corresponding to the 
separable potential V (x,t) and how this basis is easily used to express the general solution 
of the initial value problem. We then give a qualitative description of the kinds of functions 
Vo(x,t) one can obtain from this procedure. As we ultimately want to consider perturba- 
tions of Vq(x, t), we next show how to use the basis of solutions to the unperturbed problem 
to derive the coupled mode equations which trivialize the unperturbed dynamics and lay 
bare the perturbative effects. Finally, we specialize to the case of an even periodic potential 
corresponding to a two-soliton solution of the cubic nonlinear Schrodinger equation. As one 
might anticipate, the evenness (in x) of the potential leads to some simplifications in the 
spectral representation. 

2.1 Separable time-dependent potentials. 

Let us present the construction of the family of time-dependent potentials that we will 
consider in this paper, and describe their properties with respect to the linear Schrodinger 
equation. More details can be found in |Appendix A. 



Each potential we shall consider will be specified by a certain set of discrete data. Let 
TV and M be independent natural numbers. A set of discrete data T> consists of M distinct 
complex numbers Ai, . . . , Xm in the upper half-plane, and M vectors in C N . 

The discrete data V is used to build a potential function Vo(x, t) in the following way. 
Introduce the scalar expression 

a(x,t,X) = \ * M +J1 APfl(p) (M)j e~ 2l ^ Xx+x2t \ (2.1) 



and the iV-component vector expression 

b{x,t,X) = ^ P b {p) {x,t). (2.2) 



M-l 



p=0 
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In these expressions, the coefficients a,W(x,t) and b<*>(x,t) are undetermined functions of x 
and t. They will now be determined with the use of the discrete data T>. For k = 1, . . . , M, 
we insist that a(x, t, A) and b(x, t, A) satisfy the relations: 



a(x,t, A fc ) 
b{x,t,\%) 



g^b(x,t,\ k ), 
-a(x,t,\l)g^ 



(2.3) 



These equations amount to a square linear inhomogeneous system of algebraic equations for 
the coefficient functions a^(x, t) and the components of b^(x, t). We will soon illustrate this 
procedure with a concrete example. From the solution of this linear system, the potential 
function connected with the discrete data T> is given in terms of the components of M M_1 ) (x, t) 
by 

2 



N 



V (x,t) = -4J2 



n=l 



e~ i} 0M) 



(2.4) 



This function Vo(x,t) is a genuinely time-dependent potential well. Furthermore, it can be 
shown that Vq(x, t) is in the Schwartz space as a function of x, and its L 1 norm is constant in 
t. The latter follows from the fact that Vq(x, t) can be viewed as the self-consistent nonlinear 
potential for an iV-component cubic nonlinear Schrodinger equation, which conserves the 
sum of the L 2 norms of the N field components, which are proportional to the bl^~ 1 \x,t) 
for n = 1 N. 



2.2 Solution of the linear Schrodinger equation with a separable 
potential. 

Along with the potential function Vq(x, t), this construction starting from the discrete data 



T> also provides all of the solutions of the corresponding linear Schrodinger equation |p2 
These are built from the function a(x,t, A) as follows. For all real A, set 

7r IIl A_Afe lj a(x,t,X), (2.5) 

and let the functions ^^ t i(x,t), . . . , \I/b,Af(a ; , t) be the result of applying the Gram-Schmidt 
procedure (in L 2 (M) with respect to x) to the functions a(x,t, \*[), . . . ,a(x,t, \* M ) at any 
fixed value of t. Then we have p2]j : 



1. Each function ^^(x,t,X) for A G M and each function ^/^k(x,t) is a solution of the 
linear Schrodinger equation with potential Vq(x, t). The fact that the L 2 inner product 
is an invariant of the evolution shows that the functions ^^ t k{x,t) do not depend on 
the choice of the time t at which they are obtained from the Gram-Schmidt process. 

2. For any fixed t, these functions form an orthonormal basis of L 2 (IR). 

These facts show us how to solve the initial value problem for the linear Schrodinger equation 
for the potential Vo(x,t). Namely, to find the solution of 

id t f + l -d 2 J - V Q (x, t)f = 0, f(x, 0) = f (x) e L 2 (R) , (2.6) 
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one simply projects the initial data onto the basis at t = by defining: 

/(A) = (¥„(., 0, A), /„(•)> , A = 0), /„(.)) , (2.7) 

and then recovers the solution as a superposition of modes: 

M poo 

f{x, t) = J2 fk*b,k(x, t) + / f{\)* d {x, t, A) dX . (2.8) 
k=i J -°° 

If the potential Vq(x, t) is slightly perturbed, it may still be convenient to expand in this 
basis, but then the coefficients /(A) and fx, ... , /a/ will become time-dependent. 



2.3 Qualitative description of separable potentials. 

Let us describe the types of potential functions V (x,t) that can be obtained by this proce- 
dure. In the generic case when the real parts of the parameters Ax, ... , Xm are all distinct, 
these potentials have the form of a collision among M moving potential wells. That is, as 
t — > ±oo, 

M 

K,(x,t)~^\/ (fc)± (x,t), (2.9) 

k=l 

where the individual wells have the form 

V (k)± (x,t) = -4p2 sec h 2 (2p fc (x + 2a k t) - 5±) , (2.10) 

where 5^ are constants that depend on the vectors . . . , (fi M \ and where \k = &k + ipk- 
Considered in isolation, each well carries a single bound state. When the wells collide for 
finite t, the bound states can become mixed, and a state f(x,t) that is bound in a single 
well as t — >• — oo will have a component in each well as t — >• +oo. The associated scattering 



matrix can be computed exactly [21 



If some of the parameters Afc share the same real part a, then the asymptotics of the 
potential Vq(x, t) in the frame moving with velocity — 2a will no longer be stationary, but will 
be generally quasiperiodic. In particular, if all of the parameters A^ are purely imaginary, 
then the potential Vq(x, t) will generally be a quasiperiodic function of the time t. This is clear 
because taking A& = ipk with p^, real and positive ensures that the only time dependence that 
enters into the computation of V (x, t) is via the exponentials exp(±2zp|t). Such a potential 
is automatically quasiperiodic. We can further ensure that the potential function is strictly 
periodic by making the frequencies commensurate. This will be trueQ if the parameters pu 
have the form 

Pk = \/n k ^ + A, (2.11) 

where Qq is some fundamental frequency and are distinct integers. This choice ensures 
that the frequencies 

u jk = 2p 2 j - 2p\ = (rij - n k )Q (2.12) 
3 For M > 2. The potential is always periodic if M = 2 and is stationary if M = 1. 
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are all integer multiples of Qq. Only the frequency differences Uj k are important because the 
potential is given as a sum of absolute values ( [2.4|) . In fact, it can be seen from the form of 
the linear system ( |2.3| ) that 

&<*0 ( x ,t) = e 2i ^ + - +p2 ^ t G n>p ({e i ^ ht },x) , (2.13) 

where G ntP is, for each fixed x, a rational function of the exponentials exp(iu>j k t). The 
sufficiency of the relations ( |2.11| ) to guarantee time periodicity of V(x, t) with fundamental 



frequency Qo is then clear from ( |2.4|) . 

2.4 Perturbed separable potentials and coupled mode equations. 

As we have already suggested, the explicit basis of exact solutions derived in the previous 
subsection forms a natural coordinate system in which to study perturbed problems. Let 
W(x,t) be a correction to the potential energy, so that the equation becomes 

1 



i ft = --d 2 x + V (x,t) jf + W(x,t)f = H (t)f + W(t)f. (2.14) 

Here, Vo(x,t) is a separable time- dependent potential built from the discrete data V = 
{Ai, . . . , Am, {ft 1 **, ■ ■ ■ ,{ft M ^}- So we use completeness to express f(x,t) for each fixed t in 
terms of the basis of solutions of the unperturbed problem: 

/oo 
B d {t,X)^ d {x,t,X)dX. (2.15) 
■°° 

In the absence the the perturbation W(x,t), the mode amplitudes B^^ and B d (t, A), A 6 R 
are governed by the equations d t B^k = 0, d t Bd(t, A) = 0. In the presence of a perturbation 
W(x,t) coupled mode equations can be derived by projecting (|2.14j ) onto the basis elements 



^fb,k(x,t) and ^d(x, t, A). This yields the system of coupled mode equations: 

/oo 
B d {t,\)N{t,X)d\, 
■oo 

(2.16) 

/oo 
K(t,rj,X)B d (t,X)dX 1 
■oo 

for the coefficients of f(x,t), where -E>b(£) is the vector of bound state amplitudes B b!k (t), 
and where the matrix elements of the perturbation W(x, t) are explicitly given by 

M kj (t) = (y h , k (;t),W(;t)V hij (;t)) , 
N k (t,X) = (V h>k (;t),W(;t)y d (;t,X)), (2.17) 

K(t, V ,X) = (V d (;t, V ),W(;t)* d (;t,X)), 

where N k (t, A) are the components of the vector N(t, A) and M k j(t) are the elements of the 
matrix M(t). In particular, it follows that the matrix M(t) is Hermitian and the scalar kernel 
K(t,r],X) is Hermitian symmetric because W(x, t) is real. With the unperturbed problem 
exactly diagonalized in this way, this system is a useful starting point for perturbation theory. 
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2.5 Even two-soliton periodic potentials. 

In this subsection, we illustrate the procedures described above in some detail with an 
example that is important in applications and that will guide the subsequent discussion. We 
consider the case N — 1 and M = 2, and accordingly introduce the expressions 

a(x, t, A) = (A 2 + a« (x, t)X + a® ( x , t ))e~ 2i{Xx+Xh) , b(x, t, A) = (x, t) A + 6 (0) (x, t) . 

(2.18) 

Because N — 1 these are both scalar expressions, and we have at the moment four complex 
valued unknown functions, a^°\x, t), a^(x,t), b^(x,t) and b^(x,t). To find these, we 
introduce the discrete data Ai, A2, g^\ and (again, here the are complex scalars 
because N = 1). The linear equations (|2.3|) then become 



(2.19) 



(A? + a (1) (x, t) Ax + a<°> (x, f )) e -*(Ai*+A?*) = gM* (6« (x, t)X 1 + &(°> (x, t)) 
(A 2 + a (1) (x, t)X 2 + a<°> (x, t)) e - 2i ( A ^'+ A 2*) = ^(2). ( 6 (i) ( X) t ) A2 + 6 (o) ( X; t )) 

and 

b^(x,t)Xl + bW(x,t) = -g {1 \Xl 2 + a^ 1 \x,t)Xl + a {0 \x,t))e- 2l ^ x+x " t) 
b<» (x, t)X* 2 + &<°> (x, t) = -g {2) (A* 2 + a?> (x, t)X% + a<°> (x, t ))e~ 2l ^ x+x ^ 



(2.20) 



Given the discrete data V, one can solve these equations for a^°'(x, t), a^(x, t), b^°'(x, t) and 
b^'(x, t), say by Cramer's rule, and thus obtain explicit expressions in terms of exponential 
functions. 

Specializing to the case of Ai = ipi, A 2 = ip 2 (we assume without loss of generality that 
P2 > Pi), we obtain a time-periodic potential function, since the parameters A& are pure 
imaginary and then the commensurability condition is automatically satisfied for M — 2. 
Furthermore choosing = e tdl and g^ = e ld2 ensures that the potential function is even 
in x. Indeed, we then find that with s = p2 + pi and d — p 2 — p%, 

&W (x t) = 2sd Pl cosh (2p 2 x)e 2 ^^ - p 2 cosh(2 Pl x)e 2 ^+^ 

K,J d 2 cosh(2sx) +s 2 cosh(2dx) -4p 1 p 2 cos(2sdt + 6 2 -8 1 )' 

The potential function is then given by 

V (x,t) = -4\b^(x,t)\ 2 , (2.22) 

which is easily seen to be periodic in t with period L = n/(sd), and an even function of x. 
The shapes of these time-periodic potential wells are shown in Figures |T] and for Q\ = 9 2 = 
and two different choices of the parameters p\ and p 2 . From the solution of the same linear 
system, we also find 

n( ov f , nn (Pi + P2?S l S 2 + p\e-^-^) + pje^K^) - 2 Pl p 2 C 

l , J PlP2 2p lP2 cos(u;t+(9 2 -9 1 ))-(p 2 + p 2 )C + (p 1 + p 2 yS 1 S 2 ' ^ 



10 



>,= 1 .o 




Figure 1: The potential well V (x, t) for p x = l/y/2 and p 2 = 1. The phase parameters are 
9i=02 = 0. This well is time-periodic with period L = 2i\ and even in x. 



and 



a w (x,t) = i- 



(Pl - P 2 2 )piC 2 Si + {pi - p 2 1 )p 2 C 1 S 2 



'2 PlP2 cosM + (0 2 - e 1 )) - ( P i + P i)c + ( Pl + P2 ) 2 s 1 s 2 ' (2,24) 

where Sk = sinh(2pfcx), Cfc = cosh(2pfcx), and C = cosh(2(pi + p 2 )x), and where the 
frequency is uo = 2n/L = 2sd. Note that a^(x,t) is an even function of x, while a^(x,t) is 
odd. We may then write the mode function a(x, t, A) in the form 



a(x, t, A) = ((A 2 + a (0) (x, t)) cosh(-2iA:r) + Aa (1) (x, t) sinh(-2iA:r)) e~ 2ix2t + 
((A 2 + a(°) (x, t)) sinh(-2iAx) + Aa« (x, t) cosh(-2iAa;)) e" 2 ^' , 



(2.25) 



in which the first term is even in x and the second term is odd in x. Also, it is clear that 
a(—x, t, A) = a(x, t, —A). 

A particularly convenient orthonormal basis of the two-dimensional space of bound states 
is given by the formulas 



i 



V 4 (pi + pi) 
1 

V 4 (PI + P2) 



pl - Pi 

2 



-a(x,t, —ipi) + 



P2 - Pi 



-a(x, t, -ip 2 ) 



Pi Pi - P2 



a(x, t, — ipi) + 




p 2 p 2 - Pi 



a(x,t, -ip 2 ) 



(2.26) 

In this case, the even symmetry of the potential Vo(x,t) guarantees that we may choose one 
basis element to be even and the other to be odd; we are using superscripts "(e)" and "(o)" 
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= 0.25, p 



= . 7 5 



-3 



— 4 




Figure 2: The potential well Vq(x, t) for p 1 = 1/4 and p 2 = 3/4. The phase parameters are 
6i = 62 = 0. This well is time-periodic with period L = 2ir and even in x. 



to refer to even and odd functions of x respectively. These bound state solutions of the linear 
Schrodinger equation are shown in Figures || and ||. The two bound state modes are Bloch 
functions in t, with the same Floquet multiplier exp(2ip^L) = exp(2ip^L) = exp(2iplL). 
Note that, in reference to the remark made at the end of § |2.1| , the function (x, t) is 
proportional to ip = 2ib^(x,t), which is a two-soliton solution of the nonlinear Schrodinger 
equation 

' ' ' " (2.27) 



o 



Correspondingly, V (x,t) = — \ip\ 2 is the self-consistent potential. 

It will also be useful to decompose the continuum into odd and even parts. Using the 
fact that a(x, t, —A) = a(— x, t, A), define 



^\x,t,X) 



1 



^\x,t,X) 



V^WTpIW+pT) 
i 

V^WTpW^+W) 



(a(x, t, A) + a(x, t, —A)) 
(a(x, t, A) — a(x, t, —A)) 



(2.28) 



where A > 0. The solutions ^(x, t, A) and (x, t, A) are also Bloch functions in t of period 
L with Floquet multiplier exp(— 2i\ 2 L). 
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Figure 3: Equal phase contours for the even (left) and odd (right) modes superimposed on 
a density plot of the corresponding square modulus. The parameter values are p\ = 1/ \[2 
and p 2 = 1. For these values of pi and p 2 , the Floquet multiplier is equal to 1, and therefore 
these are periodic functions oft. 



These solutions of the unperturbed problem have the following inner products [22| : 

<< c) M),*S e) M,A)> = o, 

(^\.,t),¥:\;t)) = 1, (2.29) 

(^\;tA)^\;t,v)) = 5(X- V ). 

In this latter relation it is assumed that both A and r\ are positive real. Similar relations hold 
among the odd solutions, and of course everything even is orthogonal to everything odd. If 
the perturbation W(x, t) is also even in x, then this observation will allow us to treat the 
even and odd parts of the field f{x,t) in isolation to all orders in perturbation theory. 

In our subsequent analysis of the coupled mode equations for this family of periodic 
potentials, we shall assume that the perturbation W(x, t) is also an even function of x, and 
thus restrict attention to the subspace of initial conditions, f(x, 0) which are either even or 
odd in x. By the spatial symmetry of V = V + W, f(x, t) has the same parity as f(x, 0). In 
analogy with the above derivation of coupled mode equations, we can then expand f(x, t) in 
terms of (even or odd) modes of the unperturbed problem: 

poo 

f(x,t)=B^(t)¥ h a \x,t)+ / B^\t,X)^ (x,t,X)dX, (2.30) 

Jo 

where a = e if f(x, 0) is even and a = o if f(x, 0) is odd. Coupled mode equations for 
the amplitudes B^\t), B^*\t, A), A G R analogous to those derived in the absence of any 
particular symmetry can then be derived by projecting the dynamical system ( 2.14Q onto 



these even and odd basis modes. In §E| we show, by using the Floquet factorization of 
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Even Mode Odd Mode 




Figure 4: Equal phase contours for the even (left) and odd (right) modes superimposed on a 
density plot of the corresponding square modulus. The parameter values are p\ = 1/4 and 
P2 = 3/4. Here, the Floquet multiplier is not equal to 1, and the modes are not periodic in 
t, although they have Bloch form. 

the unitary evolution associated with the unperturbed dynamics, that these coupled mode 
equations can be reexpressed as the following system (c.f. the system ( |3.15|) ) which is more 
amenable to our techniques: 

POO 

id t A h {t) + 2p h A h (t) = M{t)A b {t)+ N^(t,X)A d (t,X)dX, 

Jo 

(2.31) 

id t A d (t,rj) 2rfA d {t,rj) = N^(t, V )*A h (t) + K^(t,rj,X)A d (t,X)dX, 

Jo 

where 

A h (t) = B h {t)e 2 ^ 1 , A d (t, X) = B d (t, X)e~ 2a2t , (2.32) 
and where the scalar coefficients, all periodic with period L, are 

M(t) = (M;t),W(-,t)M;t)) 
N®(t,\) = (M;t),W(;t)^ d (;t,X))e 2 ^ 2+ ^ t (2.33) 

K^(t,r],X) = (^ d (-,t,r ? ),iy(-,t)^ d (-,t,A))e 2j ( A2 -" 2 ) t . 

Remark: To avoid cumbersome formulae, we have omitted the superscripts (o) and (e), 
with the understanding that the amplitudes correspond to either one type or the other, 
depending on the parity of f(x, 0). □ 

The system (|2.31| ) may be viewed as that governing a family of oscillators: a single discrete 
oscillator whose amplitude is A h (t) coupled to a continuum of oscillators with amplitudes 
A d (t,r]), 7] e R+. 

In we shall analyze the coupled mode system (|2.31|) , and determine the detailed 
asymptotic behavior of its solutions for small W(x, t) over different timescales. 
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3 Coupled mode equations for periodic potentials. 



Consider a dynamical system of the form (|1.6|) , where both the unperturbed and the per- 
turbed potential are time-periodic with the same period L. We will encounter a concrete 
example of such a problem in §|5|. Floquet theory pj suggests the introduction of a new time- 



periodic basis, with respect to which the problem ( |2.14| ) becomes a periodic perturbation of 
autonomous Hamiltonian system. This change of basis transforms the problem at hand into 
one similar to that treated in |23J and [TBI. Similar methods are used along with resonance 



theory in a weakly nonlinear setting in [p5 | 



3.1 Floquet factorization. 

Let U(t) denote the unitary evolution operator (or propagator) of the unperturbed problem, 
so that for any L 2 (M) function f(x), f(x,t) = U{t)f(x) is the solution of the unperturbed 
problem with f(x, 0) = f(x). As a consequence of the periodicity, the evolution operator 
can be factored into two operators on L 2 (R): 



U{t) = V{t)e 



-itB 



(3.1) 



where V(t + L) = V(t), and B is independent of t. This factorization can be motivated by 
the observation that by periodicity, there is an operator Ai satisfying U(t + L) — U(t)A4, 
and that by setting t = in fact one has M. = IA(L). Since U(L) is unitary, one can find a 
self-adjoint operator B such that M. = U(L) = e~ %LB . This operator B in turn defines the 
abelian unitary group e~ lW . Now it is easy to see that V(t) = lA[t)e ltB is a unitary operator 
satisfying V(t + L) = V(t). 

Let W(t) be the operator of multiplication by the correction to the potential W(x, t), 
and set y(x,t) = V(ty f(x,t). Then, the perturbed equation ( |1.6|) becomes 



id t y-By = W{t)y, 



(3.2) 



where W(t) = V(tYW(t)V(t), a "dressing" of W(t). The form of Q is similar to the type 
of problem treated in |28j and |15|]. The "unperturbed Hamiltonian" B is time-independent, 
and self-adjoint. The perturbation \V(t) is localized, self-adjoint, and time-periodic because 
the periods of V(t) and W(t) are equal. Typically the perturbation contains frequency 
components at all overtones of the fundamental frequency, and thus the version of the theory 
described in [15] is most appropriate. The key qualitative difference between the present 
situation and that treated in [JIH] is that here the unperturbed operator can have multiple 
bound states. We will soon introduce a symmetry that removes this difficulty from the 
However, the methods of p8| and |15 can be extended to give results 



scope of this paper. However, the methods of p8f and 
on radiation damping due to the coupling of multiple discrete modes to the continuum for a 
general class of spatially localized and time-dependent perturbations [16[. 

In fact, one can simplify the problem even further by invoking the spectral theorem for 
the self-adjoint operator B. This guarantees the existence of an isomorphism V : L 2 (R) — > 
L 2 (S,c? / u) to the space of square integrable functions on some set £ with measure dfi, such 
that VB = TV where T is a real diagonal operator on L 2 (T,,dfi) (i.e. an operator of 
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multiplication by a function from £ to R). Setting z(t) = Vy(t), we find the equation 



id t z -Tz = VW(t)Vh . (3.3) 



In quantum mechanics, making the transformation from (|2.14 ) to ( |3.3| ) to facilitate the study 



of perturbations is known as going from the Schrddinger picture into the interaction picture. 

In the particular example we will analyze in detail, arising from even perturbations of the 
two-soliton even potential described at the end of §||, the operator B has a single degenerate 
eigenvalue of — 2j3b < of geometric multiplicity two. By restricting separately to even and 
odd spaces of initial conditions (possible because the potential Vq(x, t) +W(x, t) is symmetric 



in x) the problem is reduced to one which, formally, is precisely of the type studied in [|15 |. 
We may then apply the methods developed in |15[ (subject to some appropriate hypotheses) 
without modification. 

We now use our explicit knowledge developed in §|2| of the unitary propagator U (t) cor- 
responding to a time-periodic separable potential V (x, t) to find the operators V(t) and B, 
and then to diagonalize B. This effectively implements the program described above and 
casts the perturbed problem (|2.14j ) in a form (|3.3|) more suitable for analysis. We begin with 



the observation that each element of the basis of solutions of the unperturbed problem is a 
Bloch function or Floquet mode. We have: 



# d (x, t + L, A) = e- 2lX L V d (x, t, A) (3.4) 



where A is arbitrary real. Also, we have 

a(x, t + L, -ip k ) = e 2l ti L a(x, t, -ip k ) ■ (3.5) 

Note that the commensurability relations ( J2.11Q imply quite generally that the Floquet 
multipliers exp(2zp^L) are all equal. This generalizes the observation made above in the 
context of the two-soliton potentials. This means that the entire M-dimensional subspace 
of bound states consists of degenerate Floquet modes. In particular, the elements of any 
orthonormal basis {^b,k(x, t) , k = 1, . . . , M} have the same Floquet multiplier[]exp(2z/?b-^)- 
It now follows from ( |3.4| ) and ( |3.5| ) that the functions defined by 

^\x,t,X) = e 2ix2t V d (x,t,\), 

(3.6) 

are time-periodic with period L, as denoted by the superscript "(p)". 

As described in §0, the solution of the unperturbed problem with initial data f(x) is 

4 Recall that the Floquet exponents are not unique but that the Floquet multipliers are. Identification 
of the Floquet exponents with a single number /3b > amounts to a particular choice of branch of the 
logarithm. 
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expanded as 
f(x,t) = U(t)f(x) 



X)<*b,*(-, 0), f{-))* h>k (x, t) + / (* d (-, 0, A), /(0)*d(*, *, A) dA 



fc=l 

Al 



= E^£(-,0),e 2 ^V(-))*S(^0+ / M p) (-,0,A),e- 2iA2 V(-))*d P) (^^A)dA. 
fe=i 

(3.7) 

We now use the completeness relation at t — to factor tl(t) as V(t)e %tB where 



M 



e- itB m = 53<*g(-,O),e a ^*/(0>*S(^O)+ / (^ p) (- ) 0,A) ) e- 2 ^/(-))^ d P) (^0,A) C /A 

(3.s; 



fc=l 
and 



/oo 
(^ p) (,0,A)^(-))* d P) (^^A)dA. (3.9) 
■°° 

We have used, several times, the fact that at t = there is no distinction between the basis 
elements and their periodic counterparts defined by ( |3.6|) . It is easy to see that f(x,t) = 
U(t)f(x) = 7 ? (t)e~ l * B /(x) is the solution of the unperturbed initial value problem with data 
f(x) G L 2 (R) and that V(t) is periodic with period L and U(L) = e~ iLB . 
The generator of the abelian unitary group e~ is 

Bf(x) = ij t e~ m f{< 



t=Q 



M 



k=l 

M 



0), -2&/(.)>*gta 0) + / 0, A), 2A 2 /(-))^ d p) (^ 0, A) dX 



/oo 
(vI/ d (.,0,A),2A 2 /(-)) v I > d(x,0,A) C /A. 
■°° 

(3.10) 

In the last step we have dropped the superscripts "(p)" since everything is evaluated at 
t = 0. This formula for the self-adjoint operator B makes clear its spectral decomposition. 
The isomorphism V takes a function g(x) G L 2 (R) to a function A&(X) for A G R and a set 
of M numbers for k = 1, . . . , M defined by 

A d (A) = (^ p) (-,0,A),^(-)), (3.11) 

and for k — 1, ... , M, 

^b,*=<*S(-,0),p(-)>. (3-12) 
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The diagonal operator T is then simply defined by 



r 



" Ai(A) " 













2A 2 



-2/3 b 



-2/3 b 



A d (A) 



(3.13) 



It is now easy to use the definition of the unitary periodic operator V(t) and the uni- 
tary isomorphism V, along with the completeness relation to compute the dressed operator 
VW(t)V^ and thus write the perturbed problem ( |2.14| ) in the simple form ( |3.3| ). The dynam- 
ical unknowns are in the range of V, the space L 2 (H, dfi), and are given in terms of f(x,t), 
the solution of (|2.14j ), by 



A h , k = (vv(tyf(.,t)) h>k = K p )(.,t), /(.,*)>, 

A d (A) = (VP(t)t/(.,t)) d (A) = <* d p) (-,t, A), /(■,*)> 
When f(x,t) satisfies ( |2.14| ), these quantities satisfy the system 



(3.14) 



id t A h + 2{3 h A h 
id t A d (rj) - 2 v 2 A d ( v ) 



/oo 
A d (X)N {p \t,X)dX 
-oo 



(3.15) 



N^(t, v yA h + / K®(t,r],\)A d (\)d\, 

J —oo 

where is the vector of components Ab t ±, . . . , A^m, an d the time-periodic matrix elements 
are defined in terms of ( [2.171) by 

N^{t,X) 
K&(t,r),X) 



e^+^Nit, X) , 



(3.16) 



The periodicity of these matrix elements when W(x, t) is periodic with period L is also clear 
from these explicit formulas and the Bloch relations Q3.4|) and (|3.5|) for the basis of solutions; 
these imply similar ones for the matrix elements defined by ( |2.17| ). We have 

M(t + L) = M(t), 



N(t + L,X) = e- 2i ( x2+ ^ L N(t,X), (3.17) 
K(t + L, V ,X) = e 2t ^ 2 - x2 ) L K{t,T],X). 

Of course, the right-hand side of (|3.15|) is just the operator VW(t)V^ operating on the 
dynamical unknowns. Similarly, the perturbation operator W(t) operating in the space 
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-2 



can be explicitly written as 

u M 



k=l 1=1 

M 



m /"OO 

+ E/ ^| p) (i^)(*d(-,0,7 7 ),/(-))*b J ^,0)^ 

(3.18) 



fc=i 

Ai 



j " poo 

+ E/ Ar/ p) (t,A)*(M/ b ,K-,0),/(-))vI/ d (a;,0,A)rfA 

, =1 J-oo 



oo poo 



— oo ./ — oo 



For the special choice of Vo(x,t) discussed at the end of §|], evenness implies that there 
is one bound state of each parity. If the perturbation W(x,t) also has even symmetry in 
, the coupled mode system ( |3.15| ) reduces to a system of the type ( p.31| ) when the initial 



condition is restricted to either even or odd parity. It is easily checked that the unknowns 
as defined above correspond exactly to those defined in §0 for the system ( |2.31| ). 

4 Analysis of the Coupled Mode Equations 

In this section we study the structural instability of the even and odd breather modes 
introduced at the end of §||] associated with the two-soliton time-periodic even potentials. 
We first give a simple argument valid for short times that in the presence of a perturbation 
W(x,t) to the potential V (x,t), the bound state begins to decay initially. We then seek 
to capture the dynamics for longer times, primarily to show that this initial phase of decay 
does not reverse itself, but takes on a different, exponentially decaying, character. The decay 
will be first calculated formally, using asymptotic expansions and the method of multiple 
scales. Then, using the results of Kirr and Weinstein |15| , we show that at least in the odd 
case, it is possible to make statements about the decay process that are valid globally in 
time. In particular, these arguments will rigorously justify the formal results for the odd 
case, and will show that the exponential decay model is only a valid approximation until it 
becomes smaller than the dispersive part of the solution. The bound state ultimately dies 
algebraically in time, qualitatively indistinguishable from the dispersive components of the 
solution to which it is orthogonal. 



4.1 Small time analysis. The watched pot effect. 

A simple calculation carried out at the level of the coupled mode equations ( J2.31[ ) shows 
that the effect of the perturbation is to cause the bound state to decay immediately both 
forward and backward in time. More complicated calculations will be required to show that 
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the decay does not stop or reverse for longer times, although it takes on a different character. 
The approach in the small time analysis is simply to expand the solution in Taylor series: 



A h {t) = A b (0) + at + c 2 t 2 + 0{t 3 ) , 
A d (X,t) = di(A)t + 0(t 2 ), 

and use the (known) Taylor expansions of the matrix elements, in particular, 

M(t) = M(0) + M'(0)t + O(t 2 ) , 
N<P>(t,\) = N®{o,\)+O(t). 

Substituting these series into Q2.31|) , one finds: 

[id + 2p h A(0) - M(0)A b (0)] + 

2ic 2 + 2(3 h d - M(0)d - M'(0)A b (0) - 



(4.1) 



(4.2) 



d 1 (X)N^(0,X)d\ 

Jo 

id 1 ( V )-N&(0,r ] )*A h (0) = 0{t). 
Solving for d an d c 2 yields an approximation for A\,(t), valid for small t: 



0(t 2 ), (4.3) 



A h (t) = A h (0) 



i (M(0) - 2j3b) t 



- ( iM'(O) + (M(0) - 2f3 h f + r \N® (0, A) | 2 dX j t 2 + 0(t 3 



(4.4) 



It easily follows that 

\A h (t)\ 2 = \A h (0)\ 2 



I - / J / |iV (p) (0,A)| 2 dA + O(t 3 ) 
o 



(4.5) 



Note that smallness of the perturbation is not exploited in these calculations. This Taylor 
expansion shows that the initial phase of the evolution is a process of radiative decay, since 
|^b(£)| 2 < |^4b(0)| 2 for all nonzero t in some neighborhood of t = 0. The decay is symmetric 
in time. 

The fact that the decay is an order 0(t 2 ) effect is quite general^ and well-known in the 
perturbation theory of stationary Schrodinger equations. It has an interesting interpretation 

5 In the general setting, the decay is a simple consequence of the Cauchy-Schwarz inequality. One supposes 
that U(t) is the unitary propagator of the possibly time-dependent unperturbed problem: 

iU t {t)(t> =H {t)U{t)<t>° , 

for all states <f>°. One then considers the perturbed equation 

irl>t = («o(*)+W(t))V, 
by setting i)){t) — U(t)<j>(t), giving the "interaction picture" equation 
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in the quantum theory of ideal measurements, the so-called "watched pot effect". Sup- 
pose that an ideal measurement is made at some point during the evolution of the wave 
function to determine whether the state is bound, and the measurement yields a positive 
result. The probability of a positive result at time t is |v4t,(t)| 2 /|v4b(0)| 2 . The theory of 
ideal measurements says that as a consequence of the measurement disturbing the system, 
the wave function "collapses" upon a positive result to the bound state, and evolution of 
the wave function according to the Schrodinger equation continues from this "reset" bound 
state. One may then try to determine the asymptotic effect of making many such measure- 
ments in a finite time interval. In particular, we can ask about the limiting probability of 
finding the system in the bound state after each of n ideal measurements performed at times 
t n = T/n, as n — > oo. After each positive result, the wave function collapses and the exper- 
iment is restarted. The Schrodinger evolution takes place over short time intervals so it is 
appropriate to replace the probability in each interval p(t) by its short-time approximation 
p(t) = 1 — (at) 2 + 0(t 3 ). The n measurements are independent events, so the probability of 
always finding the system bound after each measurement is simply 

P n =p(T/n) n . (4.6) 

Because the "time slice" decay probability 1 — p(t) is quadratic in t, P n tends to unityQ as 
n — > oo, regardless of the value of T. So if the measurements are performed infinitely often, 
the decay of the bound state never occurs. The quantum "watched pot" never boils. 

4.2 Multiple scales analysis. 

We begin the multiple scales analysis by assuming that the correction W(x,t) = W(x,t;e) 
to the potential energy has an expansion in a small parameter, e (see, for example, equation 

W(x,t;e) = eWi{x,t) + 0(e 2 ). (4.7) 
It then follows that the coupling coefficient functions in ( |2.31| ) have formal expansions for 

SlXlctll € * 

M{t) = eM 1 (t) + e 2 M 2 (t) + 0(e 3 ), 
N^(t,X) = eiV 1 (p) (t,A) + 0(e 2 ), (4.8) 

K<P*(t,7i,\) = eK[ p) {t, V ,\) + 0{e 2 ). 

which one solves by Taylor series in t. The result is: 

m = (i- mm + l - (-iW(o) - w(o) 2 + [w (o), w(o)]) + o(t 3 )^j m . 

The probability of remaining in the unperturbed state is then found to be (using self-adjointness of both 
W(0) and Wo(0)) 

\{u(t)mM{t)m)? = iwo),#t)>i a 

= - (l|w(o)0(o)||l||0(o)||| - 1(0(0), w(o)m)\ 2 )t 2 + o(t 3 ) . 

This quantity is initially decreasing in time as a consequence of the Cauchy-Schwarz inequality. 

6 The superlinear nature of the decay probability is important. If p(t) = 1 — \at\ + 0(t 2 ), then P n tends 
to e~ |aT| < 1 instead. 
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Here, M x (t) , N[ p) (t) and K[ p \t,rj, A) correspond to the expressions for M(t), and 
in (P^) with W replaced by Wi. 

The amplitudes v4b(t; e) and Ad(t, A; e) are assumed to have asymptotic expansions of the 
form 



A h {t;e) 
A d (t,X;e) 



4 0) (T , 7\, T 2 , . . .) + eA« (T , T 1; T 2 , . . .) + e 2 ^ 2) (T , T x , T 2 , . . .) + 0(e 3 ) 
( ° W ^ ^ ,A) + e4 1) (T ,T 1 ,T 2 ,...,A) + O(e 2 ), 



A^(r ,Ti,T 2 



(4.9) 

where the Tt = £ k t are time scale variables. Such expansions of given functions Ab(t; e) and 
Ad(t, A; e) are highly nonunique. However the guiding principle of the method of multiple 
scales (see, for example, |12|]) stipulates that the dependence of the various terms on the 
"slow" times Ti, T 2 , and so on is chosen so that each term is uniformly bounded as a 
function of the "fast" time T . This procedure is quite systematic, and is supposed to keep 
the error terms in any truncation uniformly small in time intervals where Tj, is bounded for 
some ka.se tends to zero. We will see by comparison with the rigorous results that this 
formal procedure indeed works as advertised. 

One now substitutes these expansions into the system (|2.31|) and expands the time deriva- 
tive operating on the expansion coefficients in (|4.9|) according to the chain rule: 



d t = d To + ed Tl + e 2 d T2 + ... . 



(4.10) 



The coupling coefficients, being all periodic functions of t with period L independent of e, 
are taken to be explicit functions of t = To only. Substituting them into (|2.31|) along with 
the expansions ( f4.9|) and the chain rule formula (|4.10| ), and equating terms with the same 
powers of e leads to a hierarchy of equations: 



0(1) 



0(e): 



id To A h 



(0) 



idr.A^ 



+ 



+ 



2A>4 0) 

2^Af( V ) 







d Tl A { ^ + M 1 (T )A 



(0) 
b 



! / N?>(T ,\)A$>{\)d\, 
'o 



-idr^^+N^iTo^TA^ 



POO 

+ / 1^(T , 77, A)4 0) (A) dA 
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(4.11) 



(4.12) 
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0(e 2 



&r A£> + 2(3 h A[ 



-id Tl A^ -id T2 A^ 
+ M 2 (T )4 0) + M!(T )4 1} 



+ 



+ 



iV 1 (p) (r ,A)^(A)rfA 



(4.13) 



^ p) (T ,A)^(A)dA 



(0), 



I id To A%\ v ) 2 V 2 A%>( V ) 



2/1(2), 



and so on. Our initial conditions are encoded in the expansions (PD as A { ® ] (0,0,0,...) = 
A w , Af\0, 0, 0, . . . , A) = 0, and for j > 1, A^(0, 0, 0, . . .) = A$>{0, 0, 0, . . . , A) = 0. We 
now proceed to solve the hierarchy sequentially. 

Solving the equations ( |4.11| ) at order 0(1) subject to the initial conditions gives 

A^ = Ce 2 ^ To , Af\ V ) = 0, (4.14) 

where C = C(Ti,T 2 , . . .) satisfies the initial condition O(0,0,...) = A hQ but is otherwise 
undetermined at this stage. 

In the first of the two equations (|4.12|) appearing at O(e), it is natural to make the 
substitution 



.4 



(i) 



fe 



2if3 b T 



which leads to the equation 



d T J = -d Tl C-iM 1 (T )C. 



(4.15) 



(4.16) 



Integrating with the use of the initial condition f(T = 0) = 0, and keeping in mind that T\ 
and T are to be thought of as independent variables, leads to the expression for A^': 

-To 



A 



(i) 



-d Tl C-T -iC 



M x {s) ds e 2if3hTo . 



(4.17) 



We need for this correction to be bounded as a function of T so that the asymptotic expan- 
sion will be well-ordered for long times. Since Mi(s) is a periodic function of period L, this 
requirement uniquely determines d^C: 



d Tl C=-iM 1 C, Mx = 



MUs) ds. 



Thus, 



C = De~ m ^ , D = D(T 2 , . . .) , D(0, ...)=A h0 . 
Putting together what we have for the bound state amplitude at this time, 



A (0) = De -iMrT^i Ph T ) A (l) = _ iDe -iMiTi e 2ifl b T I _ dg 



T 



(4.18) 
(4.19) 

(4.20) 
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This has been the first application in our calculation of the guiding principle of the method of 
multiple scales, that dependence of expansion terms on "slow" times is chosen to ensure that 
the expansion terms are uniformly bounded with respect to the "fast" time T . Now we solve 
for the correction to the dispersive mode amplitude at this order (in fact the leading term) 
using the second of the equations (|4.12|) . Substituting the expressions from the previous 
order and using the initial conditions gives a unique expression: 

A$\ v ) = -iDe~ mTl f° N[ p \s, v y e ^ e -^ 2 (T -s) ds (421) 

Jo 

Continuing systematically with the equation (|4.13|) for the bound state amplitude correc- 
tion at order 0(e 2 ), we substitute all the expressions known thus far and observe the utility 
of the change of variables 

A m =he -m Tie 2ip h T (422) 

We find for h the simple equation 

d To h = -D • (Mi (T ) - Mi) / ° (M 1 (s) - Mi) ds - d T2 D - D • 7 (T ) , (4.23) 

Jo 

where 

poo pTq 

7 (T ) = ^M 2 (T ) + / N[ P) {T ,X) e- 2iix2+M( - T °- s) Ni P \s,\ydsd\. (4.24) 
Jo Jo 

Equation ( ^231) can be analyzed as follows. By linearity, we can express h as a sum: 
h — h\ + h,2, where 



d^h = -D ■ (M x (T ) - Mx) / ° (Mi(s) - Mi) ds . 

Jo 

d To h 2 = -d T2 D - D ■ 7 (T ) , 



(4.25) 



and where we assume the initial conditions hi(To = 0) = /i2(To = 0) = 0. Integrating the 
equation for hi exactly using the initial condition gives: 



D 



To 



hx = --^J (Mi(s)-Mi)dsj , (4.26) 

which is periodic, and in particular bounded, by periodicity of Mi(T ). 

We now want to select the dependence of D on the slow time T 2 such that h 2 , as found 
from the second of equations ( |4.25|) , is a bounded function of T . Using the initial condition 
to integrate the equation for h 2 with respect to To while holding T 2 fixed gives 

h 2 = -T d T2 D - D / 7 (s) ds . (4.27) 



Clearly, the possibility of choosing D{T 2 ) so that the expression ( [4.27| ) is bounded in To 
depends on the behavior of 7 (Tq) in the limits Tq — > ±oo. We now study 7 (Tq) for large |Tq|. 
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We can compute the s- integral in ( 4.24 ) exactly if we introduce the Fourier series for the 
periodic function n[ p \t , A): 



2-irikTo /L 



N®(T ,\) = N 1)k (X)e 

k=—oo 

Note that in terms of the Fourier coefficients of N^ p \Tq, A; e) itself (see 

iV lifc (A) = lime- 1 iV fc (A;e). 



(4.28) 

we have 

(4.29) 



Substituting the Fourier series into fl4.24|) , integrating term by term with respect to s, and 
changing variables to a = A 2 , we arrive at 



7 



(T ) = iM 2 (T )+ J2 



4 1 /a(a - a k ) 



n,k=—oo 

where the resonances a n are defined by 

a n = nn/ L - (3 h . 



(4.31) 



Note that for all terms having a k > 0, the difference of the exponentials in the integrand 
vanishes for a = a k , so there is no singularity. Moreover, the Fourier coefficients iV ljra (A) are 
by construction analytic functions of A for A in a sector including the real axis, and so the 
quantities Nl^ n {^/a) are analytic in a neighborhood of the positive real a axis. This property 
extends to the whole integrand, and we may therefore deform the integration contour away 
from the real axis in an effort to study the behavior for large |To| by a steepest descents type 
argument. 

For positive To, we deform the contour into the lower half plane. For 5 > 0, let C 5 + be 
the contour consisting of the diagonal segment from to (1 — i)5 followed by the horizontal 
ray from (1 — i)5 to — i5 + oo (see Figure |5[). We have 



C 



C 



Figure 5: The deformed integration contours Ci and Ci 



25 



7(T ) = *M 2 (T )+ if jVi^M^l! r e -2t( g - g „)T _ e 2„(n- fc )T /L1 ^ 



t 7ci 4^^ - a k ) 



y> . f N 1<n (y/t)N ltk (y/t)* c2niln _ k)To/Ld ^ 



n,k=— oo 



= 7o(ro)+7i + (To)+7 2 + (To), 

(4.32) 

so that on the new contour Ci the two integrals converge independently. The term 7 (T ) 
is periodic in T with period L and mean value iM 2 . The term 7^ (To) is also a periodic 
function of T with period L. Its mean value is given by the terms in the sum with n = k: 

' y c « 4^0- (a - a n ) 

n=—oo ^+ v \ / 

Letting 5 tend to zero does not alter the value of the integral, and then we may use the 
Plemelj-Sokhotski formula]] to evaluate the terms with a n > to find 



7o + it = * M 2 ~ iA 2 + T 2 , (4.34) 

where 

^ r \N ltn (V*)\ 2 d* ^ r\NUVa)lda 

M ~ ^ 4^-a B ) + ^ RV -io VF(<r-<r n ) ' (435) 

n=— oo v \ / n=no 

and 

r 2 = ^ |iVl ^ )|2 . (4.36) 
Finally, consider the term 7^ (To). Its time integral, calculated term by term, is 

- E / g yV^ ^ (4.37) 

_ y f NUV*)Ni*(V*)* e -H«-*n)T da 



7 This is merely the distributional identity 

(x ± iO)- 1 = P. V. x" 1 T wr£(x) . 
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The first term is independent of T (and also of 5 > 0, since the integral converges and 
the integrand is analytic). In the second term, the real part of the exponent is negative 
for T > 0, so for T large and positive, the integrand is exponentially small except in a 
small neighborhood of a = 0. This small neighborhood gives a leading contribution to the 
integrand that is O(T ), and in particular is bounded for large T > 0. 
Putting these results together, we find that for large T > 0, we have 

° 7 (s) ds = (i~M 2 - iA 2 + r 2 )T + 0(1). (4.38) 

Going back to ( 4.27|) , it is clear that choosing 

d T2 D = -(iM 2 -ik 2 + T 2 )D, (4.39) 

will lead to a solution h 2 (T ) that is uniformly bounded for all T > 0. Also note that the 
first term in 7 (To) contributes a term 

h 2M = -i I ° (M 2 {s) -M^ds, (4.40) 
Jo 

to the expression for h 2 (T Q ). We write h 2 (T ) = h 2t M(T ) + h 2 (T ). 

To find the behavior of 7 (To) and its time integral as To tends to —00, we repeat the above 
steps, this time deforming the integration contour into the upper half plane to facilitate the 
steepest descents argument. The path of integration is now Ct (see Figure [5]). The only 
difference is in the sign of T 2 ; the correct choice for a bounded solution for all T is 

d T2 D = -(iM 2 - ih 2 + sgn(To)r 2 )D . (4.41) 

Thus, the method of multiple scales gives the following approximation to the bound state 
mode amplitude: 



A h (t) = A h ^M e -^M^M,)t e iSK 2 t e -^\t\ ^ _ le ^ ( Mi(s) _ ^ ds 
'ft \ 2 ft 

/ (M 1 (s) -Mi) ds) -it 2 (M 2 (s)-W 2 )ds + e 2 h 2 (t) + 0(e 3 ] 
Jo J Jo 



(4.42) 

f (MJsi\- MTl ds\ - if 2 f ( Mn(.n) - M?l d,M + e 2 ~hM\ 4- D(f 3 ) \ . 
It is not hard to see that an asymptotically equivalent expression is just 

A h (t) = AhQe ^t e -^T 2 \t\ e ^^t e -ir M{s)ds ^ + 0(e 2^ (4 _ 43) 

This asymptotic formula is expected to be uniformly valid as e tends to zero for all \t\ < Ke~ 2 
for any constant K. 

So the behavior of the bound state amplitude under the influence of a periodic pertur- 
bation, as predicted by the multiple scale theory, is dominated by two effects, a shift in 
frequency accompanied by exponential decay. The shift in frequency is an order 0(e) effect, 
coming from M. This shift can be traced back to the influence of the perturbation directly 
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on the bound state; there is no coupling to any other modes in this term. The order 0(e 2 ) 
effects include both a further adjustment to the frequency through the quantity e 2 A 2 , the 
Lamb shift, and exponential decay through the quantity e 2 T 2 . Clearly these two numbers are 
the real and imaginary parts of the same complex frequency. Unlike the leading order phase 
shift, both of these effects are clearly due to the resonant coupling between the bound state 
and the continuum that is introduced and mediated by the periodic perturbation. Due to 
the exponential decay, the lifetime of the bound state is seen to be approximately e~ 2 /T 2 , 
which is quite long for small e. For this reason, under small perturbations of the potential 
energy the state is called metastable. 

Remark: The validity of this expansion procedure is clearly called into question if any 
of the resonances a n are very close to zero, in which case the complex frequency A 2 + iT 2 
is potentially large. The breakdown of the expansion in this case indicates the presence 
of a parametric zero energy resonance. Note, however, that in the odd case the matrix 
element N^ p \t, A) vanishes as A tends to zero, and therefore so do the corresponding Fourier 
coefficients (and in particular they vanish to leading order in e, that is, Ai jT1 (A) vanishes 
at A = for all n). This suggests that the expansion ( |4.43| ) continues to hold in the odd 
case as the parameters p\ and p 2 of the two-soliton potential are varied so as to cause a 
resonance <x„(pi,p 2 ) to pass through zero. In the even case, however, behavior possibly very 
different from that predicted by the formula ( |4.43| ) is expected if a resonance is close to 
zero. We plan to investigate this phenomenon analytically; however in this paper we will 
demonstrate the effects of parametric zero energy resonance in both the even and odd cases 
with numerical simulations. Sudden changes in the behavior of a simple model for atomic 
ionization as a parameter is smoothly varied, causing the system to pass through a zero 
energy resonance, have recently been observed and compared with experiment by Costin, 
Lebowitz and Rokhlenko 0. □ 



4.3 Rigorous analysis and infinite time results. 

The multiple-scale analysis of the preceding section leads to an asymptotic formula for the 
decaying bound state amplitude that is valid on time intervals of order e -2 . In this section, 
we will establish the validity of the asymptotic formula ( 4.43| ) in certain circumstances using 



the results of Kirr and Weinstein [[OJ. When applicable, these results also yield a detailed 
description of the solution as t — > ±oo. 

More precisely, we now study the perturbed periodic system in the form obtained by use 
of Floquet factorization of the time-periodic unperturbed Hamiltonian 7io(t): 

id t y-By = W(t)y. (4.44) 

The self-adjoint operator B : L 2 (R) — > L 2 (IR) defined in §^ can be thought of as a time- 
independent Hamiltonian, and the idea is to apply the theory of periodic (or almost periodic) 
perturbations of autonomous linear Hamiltonian systems as developed in [HJ and directly 
to the problem in this form. 

As we did in the multiple scales analysis, we will restrict attention to the special case 
of periodically perturbed even two-soliton periodic potentials. As we know, in this case the 
operator B has exactly two L 2 eigenfunctions, one an even function of x and the other an 
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odd function of x. Since L 2 is the direct sum of its two subspaces L 2 e , of even and odd 
functions, and since B leaves each subspace invariant, we may study the problem Q4.44|) 
restricted to one subspace at a time. This reduction results in an unperturbed problem with 
a single bound state, and to such problems the results described in and |28| can be 
applied without modification. 

On each subspace L 2 e ^ (M) , the operator B is explicitly given by: 

POO 

Bf(x) = (^ e ' o) (-,0),-2/3 b /(-))^ e ' o) (a ; ,0)+ / (\I>|f o) (-, 0, A), 2X 2 f(-))^'°\x, 0, A) dX , 

Jo 

(4.45) 

where the functions ^'°\x,t) and \l/ jj (x, t , A) are defined given the parameters p\ and p 2 
in §0. The hypotheses required in Jl5| of the even and odd restrictions of the operator B are 
reproduced here adapted to our application: 

(HI) The even and odd restrictions of B are densely defined on subspaces of L?x(R) and 
have self-adjoint extensions to all of L 2 eo ^(M.). 

(H2) The spectrum of B in each of L^ eo ^(IR) consists of an absolutely continuous part 
°cont (^0 = [0) °°] whh associated spectral projection p^ 6 ' - 1 and a single isolated eigen- 
value Ao = — 2/3b with corresponding normalized eigenstate ipo(x) = ^'°\x,0), so 
that 

B^o = AoVo, ll^o|| 2 = l. (4.46) 

(H3) The odd restriction of B satisfies two dispersive local decay estimates. There exist 
constants C ns and C s such that 

(a) The nonsingular local decay estimate 

\\(-)- 7/2 e- m V^f\\ 2 < C n s(t)- 3/2 ||(-) 7/2 /l| 2 , (4.47) 

holds for all / e L 2 o) (R). 

(b) The singular local decay estimate 

\\{-)- 7/2 e- m {B-2p-2iK0)- 1 V^fh < a(t)- 3/2 ||(-) 7/2 /|| 2 , (4.48) 

where k = sgn (t), holds uniformly for all p satisfying \p\ > p m m > 0, that is, the 
constant C s only depends on p m ™- 

The local decay hypotheses are established in Appendix rj . We remark here that due to a 
zero energy resonance, the decay estimates that are established in |Appendix B| for the even 
case are of the form ( [4.47] ) and ( (4.48|) but with decay rate (t)~ 1/ ' 2 (this is a sharp estimate). 
Unfortunately, this slower rate of decay precludes the direct application of the results in |]15| 
and [p8|] . On the other hand, as long as the perturbation does not create a resonance p that 
is close to zero, we can expect similar results to hold in the even case over time scales of 
length \t\ < K/e 2 , since there is no obvious difficulty with the multiple scales analysis. 
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The application of the results of |0 and [28] also requires some hypotheses to be satisfied 



by the perturbation operator W(t) and its relation to the unperturbed Hamiltonian B. The 
perturbation operator acting on L 2 , , (R) takes the form 



W(t)f(x) = (M(t)(^ o \.,0)JX-)) 

N^\t,r ] )(^°\.,0, V )J(.))dr ] j^°\x,0) 



poo 

+ L 



(4.49) 



poo 

+ / 



+ r K^(t, A, ri) (^ o) (-, 0, i,), /(•)) d V ) ¥^\x, 0, A) d\ . 



We recall that the periodic "matrix elements" in the above expression are defined in terms 
of either the odd or even modes by (|2.33|) . This operator, being periodic in t with period L, 
has a Fourier series expansion 

oo 

W(t) = e 2wikt/L W k , (4.50) 

k=~ oo 

where each operator VV k has the same form as ( [4.49j ) with the functions M(t), N( p \t, A), and 
K^(t, A, 7]) replaced by the corresponding fcth Fourier coefficients, M k , N k (X) and K k (\,r)) 
respectively. The hypotheses required in [|15|] of the perturbation adapted to this context 



arc: 



(H4) The operators W k satisfy VV-fc = W\ and 

oo 

l|Wfc||£(L2(R)) < OO, (4.51) 

fc=— oo 

where || ■ ||/:(l 2 (m)) denotes the uniform operator norm in L 2 (R). Also, 

oo 

ll|W(OIII= E (ll(-) 7/2 ^ll^«) + ll(y /2 W fe (-) 7/2 ||£(^(m))) <oo. (4.52) 

k=— oo 

(H5) The following resonance condition holds: 

oo 

r = 7T (Wo, S(B - 2a n ) W n ^ ) > , (4.53) 



n=no 



where the resonances are defined by a n = (Ao + 2nn/L)/2 = —/5b + nir/L. Here no is 
the smallest positive integer for which a no > 0. Note that since 

POO 

W n ^ = W n ¥°>°\x,0) = M n ^°\x,0) + / iV n (A)*^ c ' o) (x,0,A)rfA, (4.54) 

Jo 
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and since for a > 



POO 

(f(.),8(B-2a)f(.))= / \(VV^f)(\)\ 2 6(2\ 2 -2cr)d\ 

Jo 



|(VP c (e ' o) /)(v^)P 



4^ 

(4.55) 

The formula for T can be written as 

r = -£^^. (4.56) 

n=n,Q v 

Note that if N n (X) has an expansion in a small parameter e, with leading term linear 
in e as was assumed in the multiple scale analysis, then the leading term of the cor- 
responding expansion for T(e) is exactly e 2 r2, where 1^2 is correctly obtained by the 
multiple scale analysis and is given by ( |4.36| ). The constant T is a decay rate associated 
with the bound state of the unperturbed system. The statement that the expression 



Q4.53|) should be positive for decay to occur as a consequence of resonant coupling 
to the continuum is attributed to Fermi and is known as the "Fermi Golden Rule". 
Again, because the decay constant T is quadratic in the size of the perturbation, the 
exponential decay process is very slow for small perturbations. Thus, in the presence 
of a small perturbation W(x, t), the bound state is said to be metastable. 

(H6) There are no finite accumulation points of the resonances a n , n > no- This is satisfied 
automatically because the Fourier expansion of W(t) is that of a periodic function. 
The point here is that the results in [[HJ are more general; for example this hypoth- 
esis is satisfied by finite Fourier sums with incommensurate frequencies. Yet further 



generalizations can be found in [15 



Verifying the hypothesis (H4) would seem to require more detailed information about the 
correction to the potential energy W(x, t) than we have used thus far. We merely point out at 
this time that by elementary Cauchy-Schwarz arguments applied to the unitarily equivalent 
operators VVV„V^, one finds the estimate 



/ roo roo roo 

l|W n || £(L 2 (R)) <2J\M k \i + 2j \N n (\)\*d\ + J J \K n (X,r])\ 2 dX drj . (4.57) 

Assuming these bounds all are finite, which is really a question of the smoothness and decay 
of "snapshots" of the function W(x, t) at fixed t, we see that the first required bound in (H4) 
will be satisfied if the Fourier coefficients of the function W(x,t) in t decay faster than, say, 
1/n. This is because the other periodic contributions come from the analytic eigenfunctions, 
whose Fourier coefficients decay faster than l/n p for any p > 0. Therefore, not much beyond 
continuity in t is required of W(x, t), at least for this simpler estimate. More restrictions are 
certainly required to satisfy the second estimate of (H4). 
These hypotheses imply the following results: 



Proposition 4.1 (Theorem 2.1 of [|T5| ) Let B andW(t) satisfy the above hypotheses and 
let an odd function yo(x) be given such that (x) 7 ^ 2 yo(x) G L 2 o ^(IR). Let y(x,t) be the solution 
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of ( 4-44 ) with initial condition y(x,0) = yo(x). Then if |||W(-)||| is sufficiently small there 
exists a constant C such that 



(•)- 7/2 y(-,t)l| 2 <^(t)- a/2 H(-) Y/ ^o( 



-3/211/ v7/2, 



(4.58) 



holds for all t E 



Proposition 4.2 (Theorem 2.2 of ||15|| ) Assume the same hypotheses of B and W(t). 
Then if |||W(-)||| is sufficiently small, the solution y(x,t) of (U-44) corresponding to the 
odd initial condition yo(x) with (x) 7 ^ 2 y (x) E L? s(R) is of the form 



y{x,t) 



where 



K o) (-,0),y (-))e 2i/3b 'e- r| ' l e iAt e-^o M W' is e iri W+r 2 (t)] ^ o) (a;,0) 
+ (e- m V^y (-))(x,t)+y(x,t), 

oo 



(4.59) 



(4.60) 



and where 

• The phase correction ri(t) is uniformly bounded and 0(|||W(-)||| 2 ) , 

• The bound state amplitude error r 2 (t) is 0(|||VV(-)|||) uniformly for all \t\ < K/Y for 
all fixed K and decays for large time as 0((t)~ 3 / 2 ) , 

• The correction y(x,t) is orthogonal to the bound state: (ipo,y(-,t)) = for all t, and 



satisfies the dispersive decay estimate ||(-) 7 ^ 2 y(- 



0((*}~ 3/2 ) for large t. 



Remark: Propositions |4.1| and [4.2| would appear to say that all initial conditions decay 
exponentially and then algebraically. However, a more careful reading shows that it is 
possible for there to be a transient stage of growth, before the decay ultimately sets in. This 
is because the error terms, although small when the perturbation is small, are not uniformly 
small for all initial conditions yo(x) such that {x) 7 l 2 yo(x) ranges over the unit sphere in 
L 2 (IR). So, for each fixed perturbation W(x,t), no matter how small, it is possible to find 
an initial condition yo(x,t) that grows before it decays. This is achieved by the following 
thought experiment. Suppose the periodic perturbation W(x,t) is fixed and even in time t. 
Now pick any initial condition yo(x) so that (x) 7 ^ 2 yo(x) E L 2 (R). Proposition |4.1| guarantees 



that after a sufficiently large number N of periods, the size of the solution of ( |4.44| ) when 
measured in the weighted L 2 norm is as small as we please. Note that throughout this 
process, the solution continues to satisfy || (-) 7 ^ 2 y(-, t) || 2 < oo. So now, start again at t — 
with the new initial condition yo(x) = y(x,t = NL)*. Since the potential is real and even 
in time, integration of ( 4.44|) with this new initial condition is, up to complex conjugation, 
equivalent to integration backwards in time from the time t = NL with the initial condition 
y(x,t = NL). So we know that for this very small initial condition, the weighted L 2 norm 



32 



must first grow to an order one size at time t = NL as the decay process transiently reverses 
itself, before ultimately giving way to decay over longer times. The existence of such solutions 
does not violate the statement of Proposition fO or Proposition [O because if one keeps the 



same initial condition and then makes the perturbation smaller yet again, the connection 
with the time-reversed problem is lost for this initial condition, and decay occurs sooner. □ 
By the same arguments applied in the above discussion of the decay constant T, it follows 
that there is an alternative formula for A: 



no-l 



A= E f RV r m . ,4.61) 



n=no 



Again, if N n (X) has an expansion in a small parameter e of the form N n (X) = eAi jn (A) + 0(e 2 ) 
then the leading term of A is of the form e 2 A 2 , where A 2 as given by ( |4.35| ) was resolved 
by the multiple scales analysis. This frequency shift associated with the decay of the bound 
state is the Lamb shift. 

From these results, one recovers the true dynamics by setting f(x,t) = (V(t)y(-,t))(x,t), 
where V(t) is the periodic operator that appeared in the Floquet factorization of the prop- 
agator U{t) for the periodic unperturbed Hamiltonian Ho(t). Since 

(V(t)e 2 ^\;t))(x,t) = (V(t)e- m ^\;t))(x,t) = (U(t)¥°\.,t))(x,t) = < o) (x,t), 

(4.62) 

it follows that the time-dependent projection of f(x,t) onto the bound state Bloch function 
(x, t) is uniformly approximated by 

B h (t) = (¥°\.,t),f(;t)) ~ (< o) (-,0),/(-,0))e- r IV A V^o^)^. (4.63) 

For the system restricted to the odd part of £ 2 (1R), these theorems provide justification for 
the formal multiple scales analysis carried out above, and more. They globally describe the 
decay process for all time, where the multiple scales calculation only attempts to capture the 
dynamics over time scales of length T -1 . On the other hand, since the rate of free dispersive 
decay is not sufficient in the even case to apply this detailed theory, we must settle for the 
multiple scale expansions. 



5 Applications in Planar Waveguide Optics 

In this section, we present a physical application of the kinds of perturbed time-dependent 
Schrodinger equations we have been studying in detail. This will provide a concrete family 
of perturbations W(x, t) that we can use in subsequent numerical experiments. 



5.1 Time-dependent Schrodinger equations in waveguide optics. 

For completeness, we present here a brief derivation of the time- dependent Schrodinger 
equation as it occurs in the paraxial theory of monochromatic waveguide optics. Consider 
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Maxwell's wave equation for the electric field vector E(x, t) in a planar (x = (y, z)) dielectric 
medium with isotropic, inhomogeneous linear susceptibility x^ t) 

AE - ^E tt - V(V • E) = ^[x (1) (x, t) * E] u . (5.1) 

Here, the asterisk indicates convolution in time. A Fourier transform (denoted with the 
operator T) in t with dual variable uj (the optical frequency) leads to 

ATE - V(V • TE) + u2n2 (v> z > u ) jrg = o ) (5 .2) 

where the refractive index n is defined by n 2 (y, z, uj) = 1 + (Tx^)(y, z -, ljJ )- We now assume 
that the inhomogeneity is weak, so that gradients of {Tx ){Vi z i u ) are small. This implies 
that in the absence of any free charges, the approximate relation V • TE rs follows from 
the exact relation for the electric displacement V • TD = 0. Neglecting the divergence term 
in (|5.2|), one may then choose any unit vector e and set (TE)(y, z,u) = <f)(y, z,ui)e, which 
gives the Helmholtz or scalar wave equation for <p\ 

u 2 n 2 (y,z,u) 

<Pzz + <Pyy H — 2 <P = • (5.3) 

In the design of integrated optical devices, the inhomogeneity in the refractive index is 
a localized modulation of a "background index" no(o>). Choose a fixed length scale Lq and 
nondimensionalize by setting z/Lq = 5~ 2 Z and y/Lo = d^Y , where 5 is a dimensionless pa- 
rameter, and Y and Z are dimensionless coordinates. Setting (p(y, z) = f(Y, Z) exp(ij3Z/5 2 ), 
where (3 = L un (u)) / c is also dimensionless, one arrives at 

n 2 (YL 5-\ZL 5- 2 ,uj) 



2i(38- 2 f z + hz + 5- 2 f YY + (3 2 6~ 4 
With the definition 



Q(Y,Z;uj) = -^- 2 



n 2 (YL 5-\ZL 5- 



^-1 



Tin IUJ 



f = 0. (5.4) 



(5.5) 



we see that the formal limit of 8 I with f3 and Q(Y,Z;u) held fixed yields the paraxial 
wave equation 

iPfz + \fyy - P 2 Q(Y, Z; u)f = 0. (5.6) 

The potential function Q(Y, Z\ uj) vanishes as the refractive index approaches its background 
value n (uj), say as Y and Z go off to infinity (at least in most directions). Given a fixed 
function Q(Y, Z; u), we see that the paraxial approximation made here (5 j 0) is valid if the 
modulation in the refractive index is weak, slowly varying, and more slowly varying in the z 
direction than in the y direction. That is, a fixed function Q(Y, Z; uj) provides an asymptotic 
description of a family of physical refractive index profiles parametrized by 5 <C 1: 

n 2 (y, z, uj; 5) = n 2 (u) - 25 2 nl(uj)Q(5y/ L , 5 2 z/ L ; u) . (5.7) 

Note that these assumptions about the refractive index justify a posteriori our neglect of 
the term V(V ■ TE) in the original wave equation, because in the limit 5^0, gradients of 
n 2 (y, z, uj; 5) necessarily vanish. 
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5.2 Spectral properties of paraxial waveguides. 

In optical waveguide theory, integration (numerical or otherwise) of the linear Schrodinger 
equation (|5.6| ), also known as the beam propagation method, is one of the main tools for 
studying the optical properties of "long" planar structures like gradual fiber tapers or channel 
waveguide junctions, in which backward reflecting waves can be neglected. In this connection, 
a common problem that arises is the description of the change in behavior of a waveguiding 
structure as the optical frequency is varied in the neighborhood of some frequency ljq. If the 
structure n 2 (y , z , u) is one that admits the paraxial approximation, we can use the theory 
described above as a model. In this case, it is convenient to choose the length scale L 
so that at the frequency ojq we have [3=1. With this choice, we think of (3 = /3(u>) as a 
function of frequency satisfying f3{uo) = 1. With the function Q(Y, Z\ u) chosen consistently, 
the problem becomes one of studying the dependence of solutions of (|5.6|) on the frequency 
parameter u> near u> . With the change of variables x = Y^J f3{u) and t = Z, the equation 
( |5.6D takes the form 

ift=(~dl + V {x,t)^)f + W(x,t)f, (5.8) 

where 

V (x,t)=Q(x,t;uJo), (5.9) 
and the correction to the potential is given by 

W(x, t) = (3{u)Q{x/^{u), t; u) - Q(x, t; u ) . (5.10) 

Setting e = lo/ojq — 1, we see that W(x, t) = W(x, t; e) is uniformly small in e if Q(Y, Z\ u) 
is in the Schwartz space with respect to Y. We have the expansion 

W(x,t;e) = eoo [/3'(u ) (l - ^d^ Q(x,t]u ) + d u Q(x,P,u j\ + 0(e 2 ) , (5.11) 

uniformly in x and t. If the frequency range of interest is sufficiently small, then it is often 
a good approximation to consider the problem to be dispersionless, so that the refractive 
index n(y, z, us) is independent of u. In this paper, we will accordingly consider the function 
Q to be independent of to in which case Q(x, t; u) = Vq(x, t) for all u in the range of interest, 
and we can drop the corresponding term in ( 5.11|) . 



Suppose now that we choose to study a refractive index profile n 2 (y,z) that is even in 
y and periodic in z, such that after choosing a frequency ujq and nondimensionalizing, the 
function Vq(x, t) is one of the separable potentials described in detail at the end of §^|. Over 
length scales where the paraxial approximation is valid, this periodically modulated channel 
waveguide will actually have two "breather modes" , approximately described by the bound 
states ty^\x,t) and ^^\x,t). The effect of not being fully in the paraxial limit (that is, 5 
is small but finite) is that the modes will very slowly attenuate as they propagate forward 
due to a small coupling to backward propagating fields. This small attenuation occurs at all 
frequencies near u = Uq in a way that can be quantified [p3|] . However, the profile n 2 (y, z) is 
very special in that at the frequency u = uoq there is no coupling between the bound modes 
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and any forward propagating radiation modes. This additional coupling would indeed be 
present for "typical" ^-periodic waveguide profiles n 2 (y,z). 

In fact, the theory developed in can be applied to the perturbed problem ( |5.8|) because 
the unperturbed potential Vo(x,t) and the perturbation W(x,t;e) are both even functions 
of x that are periodic in t with the same period L. This theory shows that the additional 
attenuation due to coupling to forward propagating radiation, while completely suppressed 
at the frequency u , reemerges upon detuning the frequency slightly from uq. Suppose the 
waveguide is cleaved at z = and is illuminated at this face with a broadband source con- 
sisting of many frequencies u. After some distance all of the frequencies will have attenuated 
somewhat due to backscattering (weak non-paraxiality) . However, all frequencies except u 
will additionally decay by forward propagating radiation damping. The waveguide will there- 
fore preferentially "pass" light of the frequencyQcuo- These effects were observed numerically 
in 1,1- 

Note that from the point of view of optical waveguide theory, the periodicity of the 
index n(x, z) in z is an important feature, since it gives rise to an attenuated frequency 
response that is a symmetric function of frequency to in the neighborhood of ujq. Thus, 
attenuation occurs whether u is less than or greater than the frequency u)q of structural 
instability. By constrast, channel waveguides, where n(y, z) is independent of z, also exhibit 
frequency-dependent structural instability at cutoff frequencies where the number of bound 
states changes. However, in such waveguides the number of bound states (guided modes) 
is always an increasing function of frequency pOfl , which implies that an input beam that 
matches onto a mode at its cutoff frequency uq will attenuate for u slightly less than ujq 
but will remain bound and thus give rise to a significant transmission for u slightly greater 
than ujq. Thus, whereas channel waveguides with ^-independent refractive index profiles can 
behave as "high-pass" components, ^-periodic waveguides that at frequency ujq are modeled 
by separable potentials can behave as "band-pass" components. 

6 Numerical Simulations 

Here, we describe some numerical simulations we performed to verify the analytical predic- 
tions where we expect them to apply. We also would like to explore the behavior of the 
perturbed system in parameter regimes where we expect zero-energy resonances (see the 
remark at the end of § |4.2| ) to prevent the theory from applying in its current form. For 
concreteness, we considered periodic perturbations of two problems, each associated with 
a particular two-soliton separable periodic potential. The particular perturbation we se- 
lected was exactly the type considered in §0, namely, given a separable two-soliton periodic 
potential V (x,t), we numerically integrated the equation 

idtf + ^^fJ " (1 + e)V Q (x, t)f = 0, (6.1) 

8 Actually, the "background" attenuation due to nonparaxiality (5 ^ 0) decreases slightly with increasing 
frequency. When this effect is combined with the frequency-dependent decay calculated from the paraxial 
approximation, the preferred frequency for which the loss is minimal is detuned slightly upward |23|] by an 
amount that is 0(S 2 ). 



36 



for several small values of e. This problem differs from the type to which the theory developed 
above applies only by a rescaling of x; in particular, the time scale is unaffected. 

Let us give some details about our numerical scheme. We used a Fourier split-step method 
with a local truncation error of 0(At 3 ) f24| . The spatial domain [xl,xr] of [—80,80] in the 
"non zero-energy resonance" case and [—40, 40] in the "zero-energy resonance" case (see 
below for more details about these two cases) was discretized into 1024 points. The scheme 
splits the Hamiltonian into two parts: Ti(t) = TCi + 7i.2(t), where 

Hl = ~2iTT7) d ^ K 2 (t) = (l + e)V (x,t). (6.2) 

Let U e (t, s) denote the propagator associated with ( |6.1| ). Let U\ (t — s) and ZY|(t, s) be those 
associated with Tii and 7^2 (£)• Then, the numerical scheme approximates the true integration 
over a time step of size At as follows: 

U\t + At, t) ~ U{(At/A)U e 2 {t + 3At/4, t + At/A)U[{At/A) , (6.3) 

which has an error of order At 3 . It is easy to see that, after getting started with a quarter- 
step, and until finishing with a quarter-step, iterating this approximation to the propagator 
U e (t, s) over many steps amounts to simply alternating between U\ and IA 2 each acting over 
a half-step of length At/2. 

So, in each half-step, only one of the two parts is integrated. The half-step involving TCi is 
carried out in the Fourier transform domain where one multiplies by the explicit exponential 
of the operator. This step is thus exact in time, so that the only error appears in discretizing 
the Fourier transform and is smaller than any power of Ax if the functions to be differentiated 
are taken to be arbitrarily smooth. The half-step involving 7i 2 (t) is done exactly because we 
have explicit formulas for Vq(x, t) and it is possible to find an explicit exponential of 7^2 (£)• 
That is, we can write down a formula for the multiplication operator: 

U e 2 {t, t ) = exp + e) jT V (x, s) ds^j , (6.4) 

and use it in the code. Since the temporal gradients of Vo(x, t) can be large in some parts of 
each period and small in others, we adjusted the time step throughout the period. 

We expect the perturbation to generate radiation from the central bound region of the 
potential, and we need to remove this radiation from the problem as it moves to large \x\. 
To take care of this we used a "sponge layer" in which we effectively add a term of the form 



id ■ 



<>x i ) ! - I J 1 <>x p I - 




/, (6.5) 



to the right-hand side of ( p.lfl for a positive damping factor d and width w. These parameters 
were adjusted heuristically until it was observed, roughly speaking, that no energy was being 
artificially drawn out of the center and that no energy that was radiated outward was either 
reflected or transmitted through to the other side of the periodic domain. 

We integrated for 50 periods. In all the experiments it was arranged that the fundamental 
period was L = 2tt. We initialized the field / at t = to be a snapshot of either the odd 
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or the even mode of the unperturbed problem. Then, after integrating, we calculated the 
projection of the numerical solution onto the exact solution of the unperturbed problem, 
defining: 

B h (t) = (^°\. jt )J(.,t)}. (6.6) 

We verified the accuracy of the code by checking that for e = we had Bt,(t) = 1 to several 
digits, even in the presence of the damping in the sponge layer. Note that the function B^(t) 
is related to A^(t) by the simple relation: 

B h (t) = A h (t)e- 2i ^ . (6.7) 



6.1 Away from parametric zero-energy resonance. 

For the first experiments, we selected p\ = 1/4 and p2 = 3/4 as the parameters of the function 
Vo(x, t). It is easy to check that the period is L = 2n, and that the Floquet exponent of both 
odd and even bound states may be taken to be (3b = p\ = 1/16. Therefore the resonances 
are explicitly given by 

n 1 . . 

a n = , 6.8 

2 16 ' K ' 

none of which are equal to zero. This means that there is no parametric zero-energy res- 
onance, although in the even case there still is a zero-energy resonance corresponding to 
insufficient dispersive decay. In this case, the formula for the decay constant T makes sense 
for both odd and even parity. Furthermore, for odd parity, we have a proof that the asymp- 
totic expansion obtained previously is indeed valid. 

In Figure |j, we show plots of \og(\B h (t)\) for e = 0.04, e = 0.02, and e = 0.01 for an 
initial condition of odd parity. The numerical results are plotted with solid curves, and 
superimposed are corresponding graphs of — r|t| calculated from the analytical formula, the 
analogue of Fermi's golden rule, and shown with dotted lines. The main observation here 
is that the graphs follow the corresponding straight lines, which have slopes that scale like 
e 2 , as expected. The deviation from the straight lines appears to scale like e 2 as well, and 
to decay in time. In Figure [7|, we give corresponding plots of the argument of -Bb(^) for an 
initial condition of odd parity. In these plots, it is easy to see that the phase grows roughly 
linearly in time, with slope that is 0(e). This is the contribution to the frequency shift of 
the term M, which is indeed O(e). 

Now, we consider an even initial condition, with corresponding projection -Bb(^) onto the 
even mode of the exact solution for e = 0. Figure |8| contains plots of \og(\B h (t) |) as calculated 
from the numerical data for e = 0.04, e = 0.02, and e = 0.01 shown in solid curves. Also 
plotted are the corresponding decay curves — r|t| shown with dotted lines. Although for even 



parity there is insufficient dispersive decay for the results of (T3| to apply, the decay constant 
T (or more precisely as it is obtained in the multiple scale analysis, e 2 r 2 ) is finite because 
none of the resonances a n are zero, and we see that the multiple scale theory accurately 
predicts the rate of decay of the bound state even in this case. The plots of the phase of 
B h (t) are shown in Figure [j| Again, one sees that the rate of drift of the phase is O(e), as 
predicted by the multiple scale theory. 

The significant new feature apparently contributed by the lack of sufficient dispersive 
decay for initial conditions of even parity appears to be the quality of the deviations in 
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Figure 6: The magnitude of the projection of the solution onto the bound state. Odd mode. 
Pi = 1/4, p 2 = 3/4. Solid lines are the output from numerical simulations. Dotted lines are 
the analytical predictions. 

\Bb(t) \ from the "backbone" decay e~ r '*L Not only are they larger for fixed e than for initial 
conditions of odd parity, but they have an undulatory character that suggests the possible 
contribution of subharmonic frequencies to the dynamics. The period of the undulations 
superimposed on the decay appears to be long compared with L, the fundamental period of 
the problem, but also appears to be more or less independent of e. 



6.2 At parametric zero-energy resonance. 

As a second set of experiments, we considered a potential energy function Vo(x,t) obtained 
from the parameters p\ = l/y/2 and p 2 — 1- in this case, the period is again L = 2n, 
and the Floquet exponents of both modes are (3b = Pi = 1/2. The corresponding family of 
resonances is 

n 1 . 

». = 5 -2- < 6J > 

One of these values is equal to zero. This condition for parametric zero-energy resonance 
always goes hand-in-hand with another property of this potential, namely that the Floquet 
multipliers of both modes are equal to 1. Thus, both odd and even modes are actually 
periodic functions of t with period L. 

At a parametric zero-energy resonance, the dispersive local decay estimates fail to be 
sufficient to guarantee the applicability of the theory in jlSj , also for initial conditions of 
odd parity. However, in the odd case, the formulas for the decay constant F and the Lamb 
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Figure 7: The phase of the projection of the solution onto the bound state. Odd mode. 
Pl = 1/4, p 2 = 3/4. 



shift A are finite because there is sufficient vanishing in the numerator coming from the 
missing generalized eigenfunction at o = to cancel and overcome the weaker vanishing of 
the denominator. Plots of log(|5 b (t)|) for odd parity corresponding to e = 0.04, e = 0.02, and 
e = 0.01 are shown in Figure |TU| along with dotted lines indicating the analytical prediction 
of decay. The prediction of the theory appears to be very accurate indeed. The plots of the 



phase of B^(t) are shown in Figure |TT|. They show the frequency shift scaling like e, as we 
expect from the contribution of the term M. So it appears that for initial conditions of odd 
parity, there is little if any effect of the parametric zero-energy resonance, although the rate 
of dispersive decay is smaller here than at more generic parameter values. 

Finally, let us examine the behavior of initial conditions of even parity. For such initial 
conditions and for these parameter values, we have both a simple zero-energy resonance (as 
one has in the even case for all parameter values) and a parametric zero-energy resonance 
(as occurs only for very special parameter values). It is easy to see that both T and A are 
infinite in this case, and clearly one cannot expect the multiple scale analysis to be valid. So 
what can one expect? In Figure [12] we plot log(|-Bb(^)|) for e = 0.04, e = 0.02, and e = 0.01, 
as before. This time, rather than superimposing the straight lines — r|£|, we might try to 
compare with a "renormalized" rate of decay given by the formula for T with the term coming 
from cr = simply dropped. The straight lines calculated from the renormalized version of 
T appear dotted on the plots. We still see quite good agreement at the level of a mean drift 
of |i?b(£)|. As in the previous experiment with even parity, we see subharmonic undulations 
about this mean drift. However, a key point is that whereas previously the period of these 
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Figure 8: The magnitude of the projection of the solution onto the bound state. Even mode, 
px = 1/4, p 2 = 3/4. Solid lines are the output from numerical simulations. Dotted lines are 
the analytical predictions. 

undulations appeared to be more or less indpendent of e, in this case we note that the period 
appears to scale like e _1 . Thus, there is a "slow" dynamical process involving variations of 
the amplitude that is completely missed by the multiple scales analysis in its current form. 
We must expect that whatever rescalings are required to balance the blowing up of T in the 
vicinity of a parametric zero-energy resonance will also introduce interesting subharmonic 
dynamics on the scale of T\ = et that will reproduce the effects we are seeing numerically. 
As a final remark, the phase of B^(t), as shown in Figure [13| exhibits no particularly different 
behavior than was seen in any of the other experiments. The frequency adjustment continues 
to be dominated by the relatively large term M and is therefore order e. 

7 Conclusions 

In studying the propagation of waves in time-periodic potentials, considering the problem 
at hand to be a perturbation of a separable periodic problem is evidently as easy as, and 
in many cases more convenient than working with periodic perturbations of stationary po- 
tential problems. A particular application to the theory of periodically modulated optical 
waveguides in planar dielectric media allows one to study frequency dependent attenuation 
properties of certain optical waveguides. 

Many of the difficulties described in our paper concern the influence of zero-energy res- 
onances. These are generically not present (that is, for most separable periodic potentials, 
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Figure 9: The phase of the projection of the solution onto the bound state. Even mode. 
Pl = 1/4, p 2 = 3/4. 



as for most stationary potentials), but are always present when the potential has sufficient 
symmetry, as in the evenness considered above. Many problems would therefore vanish upon 
dropping the symmetry. From one point of view, this introduces the additional complication 
of having multiple bound states that are essentially coupled to one another by the pertur- 
bation. The study of perturbed multimode problems arises naturally in the theory of light 
propagation in optical fibers having large effective cross-sections. Some of the necessary 
modifications in the theory described in |15| are described by the same authors in ||16|| . 

Of course another point of view is to keep the symmetry, and hence the possibility of 
zero-energy resonance, and study the effect of the resonance in more detail. Our numerical 
experiments suggest that the effects of such a resonance are most dramatic when the ex- 
pressions for T and A blow up, but we also see significant effects, presumably coming simply 
from the lack of sufficient long time decay of freely dispersing waves, when these quantities 
are finite. An asymptotic perturbation theory for small e should be uniformly valid with 
respect to parameters like pi and p 2 , and we plan to investigate zero-energy resonances with 
such a goal in mind in future work. 
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Figure 10: The magnitude of the projection of the solution onto the bound state. Odd mode, 
px = l/y/2, p2 = 1. Solid lines are the output from numerical simulations. Dotted lines are 
the analytical predictions. 
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Appendix A Some theory of separable potentials. 



For completeness, we here give a self-contained description of the separable potentials for the 
linear Schrodinger equation that are connected with the soliton theory of vector nonlinear 
Schrodinger equations. However, the material is auxiliary and all needed facts are reproduced 
in the main text. The results here are not new |£| but some arguments are carried out here 
in more detail. 

Let qi(x, t), . . . , gjv(x, t) be given smooth bounded complex functions of real x and t, and 
let A be the vector space of differentiate C^^-valued functions of x and t. Let A be a 
complex parameter, and consider the two linear operators acting in A: 



Or 



Or 



-2l\ f 

-Q 



(a.i; 



-2i\E + U(q)) 



43 



Odd mode, p, = 0.707, p 2 =1.00 




50 



100 



150 

t 



200 



250 



300 



Figure 11: The phase of the projection of the solution onto the bound state. Odd mode. 
Pl = l/y/2, p 2 = 1. 



where E is a matrix whose elements are given by Eij = SaSji and 







-T 

Q 




and 



T{\q) = 



-2i\ 2 + ifq*/2 \f + id x f/2 
-\? + id x ?/2 -iq*f/2 



d t - -2iA 2 E + AU(g) + -V(q) , 



where 



V(g) 



(fq* 



d x q* -iff 



(A.2) 



(A.3) 



(A.4) 



Here q denotes the column vector of the functions qk(x,t) and denotes the N x N zero 
matrix. Along with these two operators, we consider their nullspaces, K x (\,q) C A and 
Kt(\, q) C A. For generic A, these subspaces are N + 1-dimensional, and if restricted to 
generic fixed x, t, and A span C N+1 . 

If the functions qk(x,t) are chosen just right, then the subspaces K x (X,q) and Kq-(\,q) 
may coincide for all complex A: Kx = Kq- = K. If this is the case, then the common 
nullspace will certainly be contained in the nullspace of the commutator: K C K\x,T\ ■ As is 
easily checked, the commutator [X, T] is not a differential operator, but is merely a matrix 
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Figure 12: The magnitude of the projection of the solution onto the bound state. Even mode. 
Pi = 1/V2, P2 = 1- Solid lines are the output from numerical simulations. Dotted lines are 
the "renormalized" analytical predictions. 

multiplication operator, with entries depending on x and t through the qk(x,t). Since the 
kernel of the commutator contains a subspace K of dimension N + 1 for most A, x, and t, 
this implies that the operators X and T commute. It is easily checked that the compatibility 
condition [X, T] = is equivalent to the vector nonlinear Schrodinger equation: 

id t q+ \d 2 J+ {f<t)q = . (A.5) 



It is therefore necessary that flA.5 ) be satisfied by the functions qk(x,t) if we are to have 



a basis of simultaneous nullvectors in the common nullspace K. When they exist, we can 
collect all these linearly independent column vectors into a square matrix F(x,t, A). These 
ideas admit a natural geometric interpretation in the trivial frame bundle E — > IR 2 with 
fiber GL(N + 1, C). Here, X and T are covariant derivative operators for E in the x and t 
directions, and the condition [X,T] =0 means that the curvature of the afline connection 
specified by X and T is zero. This implies the existence of parallel global sections F(x, t, A) 
of the bundle E. 

Finding a global section F(x, t, A) of E given q(x,t), (that is, a matrix of simultaneous 
solution vectors) is not always easy and for this reason, we will adopt a different point of 
view below. However, it is clear from ( A.l ) and ( |A.3| ) that, given bounded functions qk{x,t) 



satisfying (|A.5|) it is possible to develop an asymptotic expansion for F(x, t, A) in the limit 
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Figure 13: The phase of the projection of the solution onto the bound state. Even mode. 

PX = 1/V2, P 2 = 1. 



A — > 00. The expansion may be sought in the form: 

F(x, t, A) = (d N+1 + A _1 F (_1) (x, t) + A~ 2 F ( - 2) (x, t) + . . .) 



e -2i(Xx+X 2 t) qT 

I N 



(A.6) 



Here, Id denotes the D x D identity matrix, and c is a complex constant. The coefficient 
matrices F^ p \x,t) are determined recursively in terms of qi(x, t), . . . ,qiy(x,t) and the con- 
stant c by collecting powers of A in the compatible equations XF = TF = 0. There is 
some ambiguity in this expansion procedure entering as integration constants at each order. 
However, it is easy to see that 



k 



N. 



(A.7) 



regardless of the values of the integration constants. 

The implications of this compatible structure for linear Schrodinger equations that are 
of interest to us in this paper are easily stated. 



Proposition A.l Suppose that (\A.3[ ) is satisfied, and let v(x, t,X) € K be any simultaneous 
nullvector of X(X, q) and T(A, q). Let V : A — > C(x, t) be the operator of projection onto the 
first component. Define the self-consistent potential 



V (x,t) = -q{x,t) T q{x,ty 



(A.8) 
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and set f(x, t, A) = Vv(x, t, A). Then it follows that 

idtf+l%f-V (x,t)f = 0. 



(A.9) 



So, for each complex A, the function f(x,t,X) is a solution of the linear, time-dependent 
Schrodinger equation with potential (|A.8| ). Solutions corresponding to different values of A 
are linearly independent. Given functions qk{x,t) satisfying the nonlinear system ( |A.5| ), one 
can look to the common nullspace K of the linear operators X(X,q) and T(X,q) as a source 
of many solutions of the linear equation ( |A.9j ). 

Remark: Let us try to put these facts in a larger context, and incidentally give the proof of 
Proposition |A.1| . It is part of the lore of integrable systems theory that linearized evolution 
equations connected with integrable systems are solvable in terms of "squared eigenfunctions" 
coming from the auxiliary linear problems making up the Lax pair for the integrable system. 
The integrable system ( A.5 ) is the compatibility condition for the equations XY = and 
TF = 0. By a change of variables (gauge transformation) F = Gexp(— iXx — iX 2 t) the two 
equations take the more familiar form of the Lax pair for (|A.5|) [T^L Bfl: 



&.G = AG . 



<9,G = BG 



(A.10) 



where 



-*r 

r 



-iX 
-<f z'AI 



B 



(A.n; 



-iX 2 +ifq*/2 Xf + id x f/2 
-X?+id x q*/2 iX 2 I-iq*f/2 

If G a and G^ are any two simultaneous matrix solutions of the Lax pair ( [A. 10 ), and if C 
is any constant (that is, x and t independent) matrix, then by setting Q 
easily obtains the equations 



GqCG^ , one 



d x Q = [A, Q] , d t Q = [B, Q] . 



(A.12) 



Equations of this form are called Lax equations, and the elements of Q are the "squared 
eigenfunctions" . The terminology becomes accurate in the scalar case N = 1 when A and B 
are in the Lie algebra si (2). In this case the solutions G of the Lax pair can be normalized to 
be in the Lie group SL(2) and therefore have determinant one. Then, because Gp is 2x2 with 
determinant one, the elements of Q are seen to be bona fide quadratic forms in the solutions 
of the Lax pair (|A.10| ). The emphasis in the literature on the sZ(2)-specific terminology of 
"squared eigenfunctions" for the forms that satisfy the Lax equations ( |A.12|) no doubt bears 
witness to the fact that so many of the famous integrable equations (e.g. Korteweg-de Vries, 
scalar nonlinear Schrodinger, sine-Gordon) are associated with si (2) representations. 

If one introduces the splitting of a matrix into blocks: M = M D + M OD where M D 
consists of the lxl and N x N diagonal blocks of M, and M OD consists of the 1 x N and 
N x 1 off-diagonal blocks of M, and if one introduces Ao = A 
is an exercise to check that the equations ( |A.12j ) imply 



A=0 



and B 



o 



B 



A=0 



then it 



i T 
-il 



OD 



+ ^ 2 Q OD - 



-il 



[B?,Q 



ODi 



[A° D JA° D ,Q OD ] 



0. (A.13) 
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If one writes 



then one finds 



Q 



OD 



o f 
h 



(A.14) 



(A.15) 



-i(9 t /i + ^d'xh + q*q r h + q*g r q* + h<fq* 







These linear equations for g and /i are consistent with the constraint h = g* at which point 
it becomes clear that g(x, t) satisfies the linearization of the vector nonlinear Schrodinger 
equation ( |A.5| ) about a solution q(x,t). 

Consider now a particular solution q(x, t) of (|A.5|) and by adjoining a new trivial com- 



ponent qjsr + i(x,t) = 0, view it as a solution (?(x,t) of ( |A.5|) in the N + 1 component case. 
From ( |A.15|) it is easily seen that the corresponding components gN+i(x, t) and hjy + i(x,t) 
satisfy 

+ ^IdN+i + (gV)&v+i 



id t g N 



+i 



-idfh 



N+l + 



2 

-d 2 x h N+ i 







(A.16) 



+i 







where we have used the fact that <f T <f* = (fq*. Now considering the Lax pair ( A.lOj) for the 
primed potentials, it is easy to see that there exists a nontrivial column vector solution of 
both equations of the form (0jy +1 , exp(i\x +i\H)) T , and that further column vector solutions 
can then be chosen to have a vanishing last component. Taking the last column of the matrix 
solution G/3 to be this particular solution, and the first + 1 columns all to have zeros in the 
final component, we see that may be inverted in two independent blocks, and therefore 
a solution of the linearized equation is given by 

N+2 

g N +i(x, t) = Qi tN+2 (x, t, A) = exp(-i\x - iXH) ^ C kjN+2 (X)G a! i tk (x, t, A) . (A.17) 

fc=i 

Since the matrix C is arbitrary, we may view the sum above as the first component of an 



arbitrary column vector solution of the Lax pair (|A.10|) with the primed potentials <f(x, t). 
Moreover, since q' N+l (x,t) = 0, the first component of a solution of the primed Lax pair 
is also the first component of a solution of the unprimed Lax pair for the fully nontrivial 
potential q(x,t). Reversing the gauge transformation between solutions G of the unprimed 
Lax pair ([A.10Q and solutions F of XF = TF = then establishes the connection with 
Proposition |A.1| . 

So, the procedure we are using for solving the time-dependent linear Schrodinger equation 
is exactly the "squared eigenfunction" linearization of a certain A^ + 1 component nonlinear 
Schrodinger equation about a particular solution having q' N+1 (x,t) = 0. The "squared 
eigenfunctions" solving the linearized problem appear to be linear in this special case because 
for q' N+ i = the primed Lax pair becomes partly trivial, and the contribution of this trivial 
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part to the matrix Q is completely explicit (the exponential function that we remove with a 
gauge transformation). □ 

We now return to the construction of self-consistent potentials and the corresponding 
solutions of (|A.9|) . The nonlinear equation (|A.5|) is an integrable system by virtue of its 



representation as the compatibility condition of two linear problems. So there are many 
well-known ways to find functions qk{x,t) for which the corresponding linear Schrodinger 
equation can be solved. But as we are interested as much in the common nullspace of X 
and T as in the functions qk{x, t), we will now describe an effective approach to finding both 
at the same time. In this approach, the object of fundamental importance is the common 
nullspace K itself. We construct it first, with the functions qk(x,t) being chosen after the 
fact precisely so that for any basis matrix F of K, we will have ^(A, q)F = T(A, q)F = 0. 

What we know about K is that whenever it exists by virtue of the compatibility condition, 
the assumption that the functions qk(x, t) are bounded (this will be justified below) leads to 
expansions for large A of a basis for K of the form ( |A.6| ). These expansions are generally 
only asymptotic; there is no guarantee that there exists a choice of the integration constants 
such that the expansion ( |A.(j| ) converges for any A at all. However, we now suppose that 
there exist solutions qk{x,t) of the nonlinear system (|A.5| ) for which an expansion ( |A.6| ) not 



only converges in some deleted neighborhood of A = oo, but actually truncates. For such 
solutions qk(x,t), if they exist, a basis of the subspace K is given exactly by an expression 
of the form 



M-l 



F(x,t,A) = ( c\ M l N+l + ^ P F {p \x,t) 

p=0 



e 



-2i(\x+\ 2 t) q1 



I 



A' 



(A.18) 



for some positive integer M, where c is a complex constant. We have multiplied by an 
explicit factor of A to bring the sum into polynomial form. 

Since we are not considering the functions qk{x,t) to be known, we do not have the 
option of solving for the coefficient matrices F^ p \x,t) by substitution into the equations 
XF = TF = 0. We therefore must consider them to be arbitrary functions of x and t until 
we know otherwise. Without any constraints on the coefficients, we see that the different iable 
matrix functions of x, t, and A of the form (|A.18| ), for given integer values of M and N, form 



a vector space Kn,m over the complex numbers. 

The space A N>M is very large. If our claim — that appropriate solutions qk(x,t) of the 
nonlinear system (|A.5|) exist - - is not vacuous, then Ajv,m should contain many proper 



subspaces that may be identified with the common nullspace K of X(X, q) and T(A, q) for 
some q. If F(x,t, A) is of the form (|A.18|) and is a basis matrix of one of these subspaces, 



then it must be determined modulo the constant c. This means that each such subspace of 
Ajv.Af should ultimately be isomorphic to C, with the isomorphism being established via the 
constant c. 

We prepare to isolate the appropriate subspaces of A NjM by defining a set of discrete 
data. Let V denote an M-tuple of pairs (\k,(f h ^) where the A& are distinct numbers in the 
complex upper half-plane and where the are vectors in C^. From each vector (f k \ we 
build + 1 vectors in C N+l by setting 

gW = {-l,$»,g£>,...J*») T E C N +\ (A.19) 
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and for j = 1, ... , N, 

b^ = (gf\ejfeC N+ \ (A.20) 
where e} are the usual unit vectors in C . 

Definition A.l A^ A/ is the subspace of Ajv,m whose elements F(x,t,X) satisfy the homo- 
geneous linear conditions 

F(x,t,\ k )a (k) = 0, (A.21) 

/or fc = 1, . . . , M and 

F{x,t,\i)$M = 8, (A.22) 

fork = l,...,M and j = 1, . . . , N. 

It is not hard use dimension counting arguments to prove the following: 

Proposition A. 2 Let the discrete data T> be given. The set of solutions of (\A.2lQ and 
(\A.22\ ) forms a one-dimensional linear subspace o/A^m- The general solution of &A.2J\) 



and (\A.22\ ) is given by the one-parameter family of matrices (\A.ltQ , indexed by the complex 
parameter c. Thus, A^ M ~ C, with the isomorphism being established via the complex 
constant c. In particular, if c is given, then the coefficient functions F^'(x,t) are uniquely 
determined as functions of x and t, and if c is chosen to be zero, then F(x,t,X) is the zero 
matrix. 

This proposition allows us to index the elements of A^ M by the constant c which is now 
a genuine coordinate for the one- dimensional subspace A^ M . We indicate the dependence 
by writing Fx> tC (x, t, A) for the matrices in this subspace. 

This proposition is true even if homogeneous constraints less structured than ( A.21|) and 



( |A.22| ) are imposed. In order for the dimension count to come out right it is sufficient to 



choose M ■ (N + 1) arbitrary complex numbers along with corresponding constant vectors 
(the numbers need not all be distinct, as long as the vectors c belonging to each 
Afc are linearly independent) and to impose F VtC (x, t,X k )^ = Ofor all A; = 1, . . . , M- (A+l). 
The additional structure in the constraints ( |A.21| ) and ( |A.22| ) is needed for the following. 

Proposition A. 3 Let discrete data V be given, and let Fz> jC (x,t, A) G A^ M . Then 

~F^-?(x, t) = - (±F<£?{x, t)Y , k = 1, . . . , N . (A.23) 



We will have use for this symmetry property below. Its proof is simple. 
Proof of Proposition |A.3| : It is sufficient to consider the case of c = 1, since the coefficient 
matrices simply scale with c. It will be convenient to introduce the block form of the 
coefficient matrices: 

a (p) $P)T 



F (p \x,t) 



(A.24) 
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where a^(x,t) is a scalar, b^ p \x,t) and t) are A^-component vectors, and D^(x, t) is 

aniVxiV matrix. We will prove the stronger result that for all p = 0, . . . , M — 1, 

^P\x,t) = -b {p) (x,ty . (A.25) 

In this form, the equations ( A.21| ) and ( A.22 ) take the form of the system: 



M-l 



p=0 
M-l 



A ^*> + ^ A P^ D (p)^(fc>_ e -M(A fc «+Ajt) 6 ( P A = o 



p=0 
M-l 



(A.26) 



=o 

M-l 



ft + E ^ (e _a( ^^ff^ + D^) = 



p=0 



where k = 1, . . . , M. From the first and fourth equations, we can eliminate t) and 

T>®(x, t), p = 0, . . . , M - 1 in favor of the t) and c< p) 0, i). This involves introducing 

the elements of the inverse W of the Vandermonde matrix V having elements Vj k = Aj -1 , 

but it leads to two decoupled linear systems, one for the (x, t) and the other for the 
cP'{x,t). These systems are: 



M 



r=l 
M 



(A.27) 



E^ (r ~ x) = -hi, 

r=l 



where 



and 



M M 

# fcr = V kr e~ 2l ^ x+X ^ + J2 V ^Y1 W; 3 V; r e- 2i ^ x+x f^^ T g^ k > , (A.28) 

8=1 j=l 

M M 

^ = \ k t k> "EE ^8^* A* M ^ fc > . (A.29) 

8=1 j=l 

It is then clear that = — M p) * for all p. □ 

So, the emphasis has changed with respect to these matrices and their relation to the 
functions q k (x, t). Rather than determining the coefficient matrices F^'^x, t) from a given set 
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of functions q k (x, t) solving ( |A.5| ) by an asymptotic expansion procedure, we are determining 
them from the discrete data T> and a choice of the constant c. If there is to be any consistency, 
then we must still have relations between the coefficient matrices F^(x,t) of Fx>, c and the 
functions q k (x,t); in particular, we can rewrite ([A. 71) : 

q k {x,t) = -F[ M ~}\x,t) (A.30) 



c 



and use it as a definition of some functions qk(x,t) in terms of the discrete data T> and the 
constant c. Note that as long as c ^ 0, then this definition is actually independent of c 
because Fd )C is directly proportional to c. The fact that (|A.30| ) is sensible as a definition of 
the qk(x,t) is shown by: 

Proposition A. 4 Let the discrete data T> be given and let the constant c be nonzero, and 
let the functions qi(x,t), . . . ,q N (x,t) be defined (in terms of V alone) via ( \A.3Q ). This 
determines the operators X(X, q) and T(A, q). Then for any Fx> )C (x, t, A) 6 A 



v 

N,M> 



X{\, q)F ViC (\, x, t) = T(A, q)F v>c {\, x,t)=0. (A.31) 

For these q k (x,t) the columns of F-£, c (x, t, A) are generically linearly independent and there- 
fore form a basis of the common nullspace K for almost all A. 

Recall that the commutator [X, T] is, for fixed x and t, a matrix multiplication operator 
in C N+1 . Thus, the existence of the common nullspace K of X(X, q) and T(A, q) of generic 
dimension AT + 1 for these functions qk{x,t) implies the vanishing of the commutator and 
the compatibility of the two linear problems. Therefore we have 

Corollary A.l The functions qk{x,t) constructed from any set of discrete data T> satisfy 
the vector nonlinear Schrddinger equation ( \A.5] ). 

A time-dependent self-consistent potential function Vq(x, t) generated from the functions 
qk(x,t) connected with a set of discrete data V according to ( |A.£| ) will be called a separable 
potential p2| . 



Proof of Proposition [A.4j : Let F v ^ c (x, t, A) G A^ M . The proof begins with the simple 

observation that, as a consequence of the vectors and in the homogeneous relations 
( |A.21| ) and ( |A.22|) satisfied by Fx) )C (x, t, A) being independent of x and t, these relations are 
satisfied by (XF-r> tC )(x, t, A) and (TFt> iC )(x, t, A) as well. For example, with the operator X, 

{X(\q)F v , c ){x,t,\k)a {h) = X(X k ,q)(F VtC (x,t, X k )a^) = X(X k ,q)0 = 0, (A.32) 

for k — 1, ... , M, and 

(X(X,q)F ViC )(x,t,Xl)b^ = X(X* k ,q)(F v>c (x,t,Xl)b^) = X(\* k ,$0 = 0, (A.33) 

for k — 1, . . . , M and j = 1, . . . , N. The argument is unchanged if X is replaced with T. 
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Next, we examine the form of the matrix (X(X, q)Fx>, c )(x,t, A). It is straightforward to 
see that 

(X(X,q)F v , c )(x,t,X) = |a m (2z[E,F( m - 1 )]-cU) + 

M-l 

A p (<9 x F (p) + 2i[E, F (p - 1} ] - UF (p) ) + (A.34) 



p=i 



(d x F<® - UF(°)) \ exp(-22(Aa; + A 2 t)E) 



Now, as a consequence of the definition of the functions qk(x, t) in terms of the discrete data 
T> and the symmetry property guaranteed by Proposition |A.3j , the leading term vanishes 
identically, that is, 

2i[E,F^ M - 1 \x,t)} = cU{q) . (A.35) 



This, along with the fact that X(X, q)Fx>, c (x, t, A) satisfies the homogeneous conditions ( |A~3 
and ( |A.22|) means that 

X{\q)Fv, c {x,t,\) e A% jM . (A.36) 

Not only that, but for matrices in A® M the only way that the coefficient of X M can vanish 
is for the leading constant to vanish. Therefore, by the isomorphism between A^ M and C 
via the leading constant, it follows that 

X(X, q)F ViC {x, t, A) = . (A.37) 

We now consider the form of T(A, q)F-p )C (x, t, A): 

T(X,q)F v , c (x,t,X) = |a m+1 (2z[E,F( m - 1 )]-cU) + 

X M ^[E,F( M " 2 )] - UF^- 1 ) - |v) + 

Af-l , . x 

A p ( d t F^ + 2i[E, F^] - UF^" 1 ) - -VF^ J + (A.38) 

p=2 ^ ' 



X { d t F^ - UF(°) - iVFW ) + 



^F(°) - ~VF<°M | exp(-2z(Ax + A 2 t)E) . 



Once again, the definition of the functions qk(x,t) and the symmetry property of Propo- 
sition |A.3| guarantee that the coefficient of A M+1 vanishes. We shall now show that the 
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coefficient of A M vanishes as well. Begin by writing F^ M t) in the block form: 



a (M-l) b (M-l)T 
g(M-l) D (M-1) 



(A.39) 



We already know by definition of the functions q^{x, t) and Proposition |A.3| , that U M ^ T = 
(f /2i and c( M_1 ) = — &( M-1 )*. Making use of the fact that all of the terms in ( |A.34j ) vanish 
identically, we also have 

SsF^-^O, t) + 2i [E, F (M - 2) (s, t)] - U(g)F (M - 1} (a;, t) = . (A.40) 

This implies that for the coefficient of A M in ( |A.38| ) to vanish, it will be enough to show that 

cV(q) = 2id x F {M ~ 1) (x, t) . (A.41) 



From ( A.4|) , it is clear ( |A.41| ) is satisfied in the off-diagonal blocks. To show that the diagonal 
blocks also vanish, we write out the diagonal blocks of ( |A.40| ): 



d x a( M -^ 




T 

a.D^^ 1 ) 







-? 



T 







(A.42) 



Eliminating the derivatives of a^" 1 -* and D^ M_1 ^ between this equation and the diago- 
nal blocks of ( |A.41| ) and comparing with the definition (AA) of V(g), we finally see that 
( |A.41| ) is satisfied identically. By similar arguments as we used above, it follows that 
T(A, g)Fx> jC (x, t, A) is also the zero element of A^ M . This ends the proof of the proposi- 
tion. □ 

We now return to the problem of interest, namely the algebraic construction of separable 
time-dependent potentials for the linear Schrodinger equation and of a large number of exact 
solutions to this linear equation. From the construction of the subspace A^ M we can extract 
a simpler construction of the quantities of immediate interest, and cast the whole procedure 
in the form of an algorithm. The key observation is that it is sufficient to build from given 
discrete data V only the first row of a matrix Fx), c (x,t, A) in the space A^ A/ . This gives us 
both the functions qk(x,t) via the first row of the coefficient matrix F^ M ~^ (x,t) from which 
we find the potential Vq(x, t) and also the image of the projection operator V that consists 
of solutions of the linear Schrodinger equation with this potential. 

So we consider the first row of Fx> c (x, t, A) and impose the homogeneous linear constraints 
( [Oil) and (TA32D . Introducing 



a(x, t, A) = Fii (x, t, A) = A 



M 



M-l 
p=0 



ip). 



X, t) 



-2i(\x+X 2 t) 



and 



M-l 



b(x, t, A) = (F 12 (x, t, A), . . . , F hN+1 (x, t, A)) T = XP b (p) (x, t) 

p=0 



(A.43) 



(A.44) 
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the relations ( A.21|) and ( |A.22|) take the simple form: 

a(x,t,X k ) = (fiWb(x, t, A fc ) , (A 451 

b(x,t,X* k ) = -a(x,t,\* k )gW, 

where k — 1, . . . , M. Note that without loss of generality, we are taking c = 1. Written out 
in its entirety, this is a square linear system for the M ■ (N + 1) unknowns, t), and 

the N elements of b^(x, t) for p = 0, . . . , M — 1. The matrix of this system, and the right- 
hand side, are explicit functions of x and t through the exponential functions contributed by 
a(x, t, Afc) and a(x, t, \* k ). 

From the solution of this linear system, one computes the potential function as: 

V (x,t) = -Aj2\b { n M - l) (x,t) . (A.46) 

n=l 

Then, we see that a(x, t, A) and all the elements of b(x, t, A) are solutions of the linear 
equation 

id t f + ~d?J-V o (x,t)f = 0, (A.47) 

for fixed but arbitrary A. Being polynomial in A, each element of b(x, t, A) sweeps out an 
M-dimensional space of solutions as A varies. The solutions contained in a(x,t, A) are more 
interesting because the presence of the exponential means that all of these solutions for real 
A are linearly independent. This immediately gives an infinite-dimensional space of solutions 
to the linear Schrodinger equation. 

In fact, the function a(x,t,X) contains an L 2 (R) basis of solutions of the Schrodinger 
equation as the parameter A is varied p2| . In particular, the set 

{a(x, t, A*), . . . , a(x, t, X* M ), a(x, t, A), A real} (A. 48) 

considered as functions of x for fixed t, is complete. For real A, set 

(M \ - 1 / 2 

7r IIl A - Afc | 2 ) a(x,t,X). (A.49) 

The subscript "d" indicates solutions that superpose to form dispersive waves. For A and rj 
real we then have the inner products 

(* d (;t,X),* d (;t, V ))=5(X-ri), (A.50) 

and for k — 1, ... , M, 

(a(.,t,Aj;),* d (.,t,/i)) = 0. (A.51) 

Also, dim span {a(x, t, X* k ), k = 1,...,M} = M as functions of x for fixed t. So, let 
{^b,A:(^,^)} be any basis of span {a(x, t, X* k ), k = 1,...,M} that is orthonormal with re- 
spect to the inner product (say obtained by the Gram-Schmidt procedure), so that 

(*bj(-,*),*b,k(-,*)) = <Sj*- (A.52) 
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The subscript "b" indicates solutions that are bound and have finite energy. Note that this 
basis remains orthonormal because the time evolution of these functions under ( |A.47|) is 
unitary. 

The completeness relation is generalized to L 2 (IR) from that proved in [^2| as: 

Proposition A. 5 Let discrete data T> be given and let t EM. be fixed. For all f(x) G L 2 (M.), 
we have the expansion 

/ d (A, t)^ d (x, t, \)d\ + Y, /b,fc(*)*b^(x, t) , (A.53) 
■°° k=i 

where the expansion coefficients are given by: 

/ d (A, t) = (M; t, A), /(•)) , U(t) = (*b,*(-, t), /(•)> • (A.54) 

The orthogonality relations for the functions ^^(xjt) and ^d{%, t, A) are implied by this 
result. Note that if / = f(x, t) satisfies ( |A.47| ) then the expansion coefficients are independent 
of t and can be constructed from the initial data f(x,0). Thus one solves the initial value 
problem for ([OTP in L 2 (R). 



Appendix B Dispersive Local Decay Estimates 

Here, we establish several important properties of the unitary group e~ ttB . We will consider 
even perturbations of the even two-soliton periodic potentials, so we will work in either the 
even or odd subspace of L 2 (M). For a given function /(•) in L 2 (M), the operator 'P^ e '°^ is 
defined as the spectral projection onto the continuous part of the spectrum of B: 

POO 

(Vt o) f)(x)= / (* d e ' o) (-,0,A),/(.))* d e ' o) (^0,A)dA. (B.l) 
Jo 

As we will now see, the main difference between the even and odd cases is in the rate of 
dispersive decay and the difference can be directly traced to the behavior of the dispersive 
eigenfunction ty^'°\x } 0, A) in the vicinity of A = 0. It is easy to see from the explicit formulas 
that the eigenfunctions are continuous in A at A = 0, and that the odd mode vanishes there: 

^ o) (x,0,A = 0) = 0, (B.2) 

while the even mode does not vanish, but is simply finite at A = 0. We say that the existence 
of a nontrivial eigenfunction at A = 0, as in the even case, indicates a zero-energy resonance of 
the system. The ubiquitous effect of a zero-energy resonance is to alter the rate of dispersive 
decay in the system. However, more dramatic effects can appear if under the influence of a 
perturbation, the zero-energy resonance is directly excited. This latter situation we refer to 
as a parametric zero- energy resonance. A system with a zero-energy resonance is "primed" 
to feel the effects of a parametric zero-energy resonance in the presence of an appropriate 
perturbation. 
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B.l Nonsingular local decay. 

First, we will prove the nonsingular local decay estimate for the unitary group e~ ttB '. 
Proposition B.l Fix a > 5/2. There exist constants > such that 

1/2 n<r/(-)ii 2 , (b.3) 
3/2 n<r/(-)ii 2 , (b.4) 

for all f G L 2 (R) for which the right hand side makes sense. 

The proof is based on a sequence of intermediate results. First, from the simple chain of 
estimates: 

||(r . e -^(e,o) /( . )||2 < ||e -^ p (e,o) /( . )||2 = r (e,o) /( . ) || 2 < y { . )h < || <.)*/(.) || 2 , ( B .5) 

we have 

Lemma B.l For all a > 0, we have the simple estimate 

||(.)-e-^ e '°)/(-)||2<||(r/(-)ll2, (B.6) 

for all f G L 2 (R) for which the right-hand side make sense. 

We now want to refine the above uniform estimate to include a multiplicative factor of 
HO^/Olh that decays in \t\. To this end, we fix t ^ and observe that by the definition of 
the operator £>, 

/oo 
( y yf(y)h(x,y;t)dy, (B.7) 
-oo 

where 

h(x,y;t) = (x)-°(y)-° / ^'°\y, 0, A)*^ e ' o) (x, 0, A)e" 2iA2 'rfA . (B.8) 

Jo 

We note here that the integral in the definition of h(x, y; t) is improper; the integrand is not 
absolutely integrable, and the integral from zero to infinity should be interpreted as the limit 
of the integral from zero to R as R j oo. This limit exists as long as t ^ 0, and consequently 
the function h(x, y; t) is well-defined for t ^ 0. The trouble with the function h(x, y; t) at 
t — is not our concern here because we already have a uniform estimate that holds for all 
t, and in particular for t near zero. Thus we will be thinking of t as being large in what 
follows. 

In any case, by Cauchy-Schwarz, we have 

\{x)-°(e-« B V^f)(x)\ < \\h(x,.;t)\\ 2 || (.)*/(•) || 2 . (B.9) 

It follows that 

||(.)-(e-^°)/)(.)|| 2 <||/i(-,-;0ll2 ||(r/(-)l| 2 , (B.10) 



and 

WW (e~ m V^f) (-)h<L^(t)- 
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an estimate that involves the Hilbert-Schmidt norm of the kernel h(x, y; t) for each fixed 
t. The rest of our work will be to show h(x,y;t) is in L 2 (R 2 ) for each fixed t, with norm 
decaying in \t\. 

First note that from the explicit formulas: 

, r (Q) fx Q x) _ 2Afl(1) (^ °) CQS (2Ax) - 2i(A 2 + qWQr, 0)) sm(2\x) 

v / 2vr(A 2 + p?)(A 2 + p|) 

, T ,(e) fT a \\ _ 2(A 2 + qW(x, 0)) cos(2A:r) - 2jAog)(x, 0) sm(2\x) 
d( " )_ V^(A' + p?)(A' + p2) ' (R12) 

where (x, t) and (x, t) are bounded analytic functions of x, we obtain 

Lemma B.2 Let the parameters p\ and p 2 be fixed. The function defined by 

q(\) = ¥^\x,0,\)¥^\y,0,\y (B.13) 

is in C fc (R + ) for all k > 0. In particular all derivatives with respect to A are uniformly 
bounded functions of X. The norms Wq^i-) ||oo are homogeneous polynomials in \x\ and \y\ 
of degree k, with nonnegative coefficients that depend only on p\ and p 2 . Also, in the odd 
case, we have q(X) = 0(X 2 ) for X near zero, while in the even case q(X) = 0(1). 

In showing that h(x,y;t) is L 2 (R 2 ) with norm decaying in t, we will find that the main 
contribution for large t comes from the part of the integral near A = 0. To see this, we first 
separate the contributions near and away from zero. Let <7a(A) be a nonnegative "bump 
function", infinitely differentiable for real A > 0, identically equal to 1 for < A < A/2 and 
identically equal to zero for A > 3A/2. Let (?a(A) = 1 — <?a(A). Then 

h(x, y; t) = h A (x, y; t) + h A (x, y; t) , (B.14) 

where 

r-3A/2 



and 



h A (x,y;t) = (x)-°(y)-° / ¥^\y, 0, X)*^°\x, 0, X)g A {\)e-^ H dX , (B.15) 

Jo 

POO 

h A (x,y;t) = (x)-°(y)-° / ^(y, 0, A)^ o) (a;, 0, A)^ A (A)e^ A2 ^A . (B.16) 

J A/2 

First, we will show that away from A = 0, we can get arbitrary decay in time. 

Lemma B.3 Fix L > 0. For some k > 2, suppose that /(A) is in C n ([L, oo]) for all 
n = 0,l,...,k. Suppose that f(L) = f'(L) = . . . = / (fc_1) (L) = and that the limit 

lim / f{X)e- 2tXH dX (B.17) 

R]oo J L 
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exists for t ^ 0. Then 



/(A)e 

where the operator A is defined by 



lim 

R]oo 



R 



2ixH dX 



< 



1 



W)(A) 



L-4 k -\t\*7>L 

d ff(xy\ 



sup|A 2 (^l fc /)(A)| , 



A / 



;b.is) 



(B.19) 



Proof: Integrating by parts k times, 



lim 



R 



/(A)e 



-2i\ 2 t 



d\ 



lim 

_Rtoo 



k—l / . \ 



ra=0 



n+1 



l) n A _1 e~ *(.A n /)(A) 



\=L 



(B.20) 



+ 



4f 



(^/)(A)e 



-2i\ 2 t 



dX 



The boundary terms at A = L vanish identically, and those at A = R tend to zero as R j oo. 
These facts prove the existence of the limit of the integral in the second line, and we find 



lim 

R]oo 



/(A)e 



-2iX 2 t 



dX 



1 



A k ■ \t\ k 



R 



lim 

R1<x> J L 
-R 



(A k f)(X)e 



-2i\ 2 t 



dX 



< 



< 



4> 



l w V Z f L |AVV)(A)|-^ (B.21) 



L-4 k - \t\ 



sup|A 2 (^l fe /)(A)| • 



\>L 



The bound is finite for k > 2. □ 

We can now apply this result to estimate hA(x,y;t). 

Lemma B.4 Fix an integer k > 2, and let a > k + 1/2. Then, the function h&(x, y; t) is in 
L 2 (M. 2 ) as a function of x andy, with norm decaying as \t\~ k . 

Proof: We apply the above lemma with L = A/2and/(A) = ^'°\x,0, X)^'°\y,0, X)*g A (X). 
This gives the pointwise estimate 



\h A (x,y;t)\ < 



i(x)- a (y)- a SUD 

A-4*.|t|* A S > U A P / 2 



x 2 (A k ¥t°\^UrM' 0, (y,v,r~gA(-))(X) 



( e .°) / 



(B.22) 



The operator A k acting on the right-hand side makes the supremum bound a polynomial in 
\x\ and \y\ of degree k. Therefore for h A (x, y; t) to lie in L 2 (R 2 ) as a function of x and y, it is 
sufficient to take a > k + 1/2. The claimed time decay of the L 2 norm is then obvious. Note 
that each derivative of #a(A) contributes a factor of order O^A^ 1 ), so the overall bound on 
the L 2 norm of hA(x,y;t) scales like A~( fe+1 ). □ 

Now, we move on to consider the part of h(x, y; t) contributed by the neighborhood of 
A = 0. We again need some technical lemmas. 
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Lemma B.5 For all /j, EM., 



Proof: First, note that 



J o 



< V3. 



< / \d(\ = \n\ 



This estimate is useful for bounded ji. Suppose ji > M > 0. Then, 



Jo 



<M + 



J M 



^ d( 



(B.23) 



(B.24) 



(B.25) 



Changing variables to r = ( 2 and integrating by parts, one finds 



M 



+ 



4/x 



AM 8 



e -2ir T -3/2 dT 



M 2 



< 



AM 



Therefore, for /x > M > 0, we have the estimate 



-2iC 2 



d( 



< M + 



AM 



;b.26) 



(B.27) 



The right hand side takes its smallest value, a/3, for M min = a/3/2. Since for < fi < M min , 
we have 



/ 

Jo 



-2i( 2 



< \fi\ < M min < 2M min = V3, 



[B.28) 



the lemma is established uniformly for all positive /i. By symmetry, the same estimate holds 
for /j, < 0. □ 

Lemma B.6 Fix L > and suppose /(A) is £mce continuously differentiable, with /(0) = 
/'(0) = 0, and /(L) = f(L) = 0. Then 



f\X)e- 2lXH d\ 



< 77TR7? SU P 
4|t| 3 / 2 <A<L 



9 2 



dX 2 \ X 



(B.29) 



Proof: Integrating by parts using the boundary conditions (evaluations at the lower bound- 
ary of A = are interpreted in the sense of the limit A I 0, that is, from above), we have 



Jo 



-2i\ 2 t 



dX = — 
At 



A d\ 



')*-ifs(f)^<«> 



Write 



-2iA 2 t _ ® I „-2ia 2 t 



dX 



J 

Jo 



da , 



;b.si) 



and integrate by parts again making use of the boundary conditions (with the same caveat 
as above), to find 

J! dX = h * dX ■ < B - 32 ' 
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With a change of variables to £ = |t| 1//2 a, this becomes 



/(A)e' 



4t • |i|V2 



L <9 2 /7(A) \ /• |4|1/2A 



<9A 2 



X 



(B.33) 



Estimating the A integral in the obvious way using the uniform bound of the £ integral by 
a/3 establishes the claimed estimate. □ 

Without the vanishing boundary conditions at A = 0, one gets less decay in time. 

Lemma B.7 Let f(X) be absolutely continuous < A < L, so that /'(A) G L l ([0,L]). Then 



/(A)e 



-2i\ 2 t 



dX 



o 



< |/(0)| + 2 / |/'(A)| dX 



V3_ 

Itl 1 / 2 



an order 0(\t\ 1 ^ 2 ) bound. 

Proof: Separate off the slow decay by writing 



f(X)e^ xH dX = f(0) 



o 



-2i\ 2 t 







dX+ / (/(A)-/(0))e 



-2i\ 2 t 



dX = I A + I 



B 



The first integral is easily transformed: 

Ia = /(0) [ L 
Jo 

and therefore easily uniformly estimated 



e~ 2tXH dX 



t l/2 



Lt 1 / 2 



-2i( 2 



d(. 



\Ia\ < 



V3|/(0)| 
Itl 1 / 2 



In the second integral, one integrates by parts to find 



IB 



f\X) / e'^dpdX 



Therefore, 



l B\ 



< sup 

0<A<L 



-2iii 2 t 



dfi 



\f'(X)\dX 



(B.34) 



(B.35) 



(B.36) 



(B.37) 



(B.38) 



< 



< 



-2iii 2 t 



dfi 



+ sup 

0<A<L 



-2in 2 t 



dp, 



\f(X)\dX 



2VS 

W 2 



\f'(X)\dX. 



(B.39) 



Combining the estimates for I a and Ib establishes the claimed result. □ 

We now want to use these results to estimate h^(x,y;t). To do this, we want to ap- 
ply Lemma [B^ or Lemma [K7| with /(A) = ¥°'°\x, 0, A)*J '°\y, 0, X)*g A (X). Now, from 
Lemma |B.2| , it is clear that the hypotheses of Lemma [B.6| concerning the behavior of / at 
A = will only be satisfied in the odd case. Here, we obtain the following. 
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Lemma B.8 Consider the odd case, and let o > 5/2. Then h,A(x,y;t) is in L 2 (R 2 ) as a 
function of x and y with norm decaying like |t|~ 3//2 . 



Proof: We have the pointwise estimate 

3y5A(s)-"(y)-" 
|/iA(a;,2/;t)| < sup 

OR ' 0<A<3A/2 



«9 2 /^ o) (o;,O l A)*i oi (y,0 > A)^A(A) 



<9A 2 V A 



. (B.40) 



From Lemma |B.2| we have that the right-hand side is a quadratic polynomial in |x| and 
\y\. Therefore for hj±{x,y\t) to be in L 2 (M 2 ) as a function of x and y it is sufficient to take 
a > 5/2. The time decay of the L 2 norm is then obvious. Note that each derivative of 5a (A) 
contributes a factor that is 0(A _1 ) so the bound on the L 2 norm scales like A -1 . □ 

In the even case, we have a zero-energy resonance, and this means that the integrand 
near A = is not small enough to allow decay as rapid as in the odd case. In this case, we 
can only apply Lemma B/7 to find the following. 

Lemma B.9 Consider the even case, and let a > 3/2. Then h&(x,y;t) is in L 2 (R 2 ) as a 
function of x andy, with norm decaying like \t\~ 1 ^ 2 . 

Proof: Using /(A) = ^ e) (x, 0, \)^\y, 0, A)*# A (A) and L = 3 A/2 in Lemma [B/7|, we have 
the pointwise estimate 



3 A/2 



\h A (x,y;t)\< I 1/(0)1 + 2 ^ |/'(A)|dAj . (B.41) 

Since the derivative with respect to A results in at most linear growth in x and y, taking 
a > 3/2 is sufficient to ensure that h&(x,y;t) is in L 2 {R 2 ) as a function of x and y. Clearly, for 
large t, the L 2 norm is 0(|t| -1 / 2 ). Note that the estimate is also 0(A _1 ) due to differentiation 
of the bump function #a(A). □ 

In both odd and even cases, the contribution of h^(x,y;t) to the L 2 norm of h(x,y;t) 
dominates for large time that of h^(x, y; t), for which we had arbitrary decay. According to 
Lemma [B.4| , for a > 5/2 this latter decay is at least as fast as |t|~ 2 . These results imply the 
following. 

Lemma B.10 Fix a > 5/2. Then, for t sufficiently large, we have the estimates: 

||(.)--( e -»^W/)(0|| 2 < j^llOVOIb, (B.42) 

and 

\\(r a (e- UB V^f)(-)\\ 2 <^J(-rf(-)h, (B.43) 



|t| 3 / 2 ' 



where K^ e > and are some positive constants. 



This result, taken together with the elementary time-independent bound established in 
Lemma [B.l| completes the proof of Proposition [B.l| . 
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B.2 Singular local decay. 

Now we prove the singular local decay estimate for the unitary group e~ ltB . 

Proposition B.2 Let \p\ > /i rmn > 0. Fix a > 7/2. Let t = nr with r > and k 
Then, there exist constants M* 6 ' - 1 > ; such that 



±1. 



oJ.U 



< M^{r)- 1 / 2 \\{-Yf{- 



12: 



(B.44) 



and 



■)- a Bm((B - 2/i - 2m5)- 1 e-^P c (°)/)(- 

oJ.U 



< M (o) (r) -3/2 ||( . r/( . )||2> 

T/ie constants M^ e '°> depend only on /i mm; so i/ie bounds are uniform for large 



(B.45) 



The proof of this proposition begins with a representation similar to QB.7|) , 

/oo 
(y)*f(y)k{x,y;t)dy. 
-oo 

where 



(B.46) 



fc(x,3/;t) = (a;)- <J (y)- <T ]im 

S W JO 



(e,o). 



x,0,A)^ e ' o) (y,0,A)* _ 2iX * 



2A 2 - 2/i - 2m5 



(B.47) 



Perhaps despite appearances, the kernel k(x,y;t) is somewhat more amenable to analysis 
than the kernel h(x, y; t) that appeared in the nonsingular case. This is because for each 
finite 5 the integrand is absolutely integrable as a consequence of the uniform boundedness in 
A of (x, 0, A) as guaranteed by Lemma p.2| and the large A behavior of the denominator. 

Again, the goal is to show that the kernel k(x, y; t) is in L 2 (IR 2 ) as a function of x and y, 
with norm that is decaying in time, although in this case we will only obtain the decay for 
t of a particular sign. First, we show that the L 2 norm exists and is finite near t = 0. 

Lemma B.ll Fix a > 3/2, [i with > // mm > ; and t with \t\ < T. Then there exist 
constants C^' ^ > depending on jj, min and T such that 



■;t) 



:ba8) 



Since the bounds only depend on \i via /i m i n , they are uniform for large \fj,\. 

Proof: Begin by setting /(A) = ^'°\x, 0, X)^'°\y, 0, \)*e~ 2tX2t . First we consider fi < 
— /-train < 0, in which case we have 



k(x,y;t) 



(x)-°(y)- 



A 2 -fi 



(B.49) 
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since there is no singularity for [i < 0. We immediately get the pointwise estimate 

(x)- a (y)- a f°° dX 



\k(x,y;t)\ < MZM^ sup | /(A) | f 

1 A>0 Jo 

ir(x)-°(y) 



o A 2 



. — sup|/(A)| (B.50) 

< ^MZMZ sup|/(A)| . 

4 V — A*min A>0 

Now we consider \i > fj, min > 0. Pick some positive G less than y//i m i n . Then 

k(x,y;t) = {X) \ l ns + I s ) , (B.51) 

where 

>Vfi-G f( X )d\ f°° f(X)d\ 



Jo A 2 - n J 



L = lim 



^Jl+G A 2 - \l ' 

/(A)dA 



;B.52) 



<5|0 J^.g (A + vV + *«5)(A - V/" + m5) ' 



where the principal branch of the square root is understood, so that the square root is nearly 
a positive number for S small. It is easy to find 

d\ f°° dX 



\U < sup|/(A)| / — — + / 
A>0 \Jo /i-A J 



Vm+g A 2 - n 

< SU ^> ol/(A)l f arctanh f + arctanh ( ) 

/ 1 / ( y/fi \ fy/U-G 

< sup f(A sup — — arctanh — f — — + arctanh 

- a>o v VmIAv/A \^+gj V 



;b.53) 



For the singular part, we find 



lim f{y/P) /^ +G + r^ +G ( /(A) _ /(VjjU rfA 

<Ho 01 + V/W + wc5 A - V A + 2^1 J X-y/Ji 

i™fWfi + /^ +G f /(A) _ ^A 



(B.54) 



64 



Therefore, 



< 



< 



jjVg)l 

2^ 
2^/2 



+ 2G sup 

|A-v^|<G 

+ 2G sup 

|A-vp|<g 



/(a) 



dA 



/(A) 



;b.55) 



A + ^ 



Again, the bounds are uniform in \i for large /x. 

Now, we simply note that the pointwise bounds for I ns and I s are themselves bounded 
by functions of x, y, and t that grow linearly at worst, as a consequence of differentiation 
of /(A) with respect to A (c.f. Lemma [B.2| for the growth in x and y, while the growth in t 
comes from the factor e~ 2tX *). Thus, to have k(x, y; t) G L 2 (R 2 ) as a function of x and y, it 
is sufficient to take o > 3/2, and then the norm will be bounded by a linear function of \t\, 
and therefore uniformly for \t\ < T. The bound is also uniform in /i for \fi\ > /i min > 0. □ 

We now note that proving the decay for large r = \t\ for // < — /U m i n < amounts to 
recalling the nonsingular local decay estimate. The integral is not really singular: 



lim 

<5|0 



-¥^\xAmr'(y,o,\) 



,(e,o), 



-2i\ 2 t 



X- 



ik5 



^ o) MX)^°>(y,0,X) 



(e,o). 



Jo A 2 - /x 

(B.56) 

Using the same arguments as used to prove the nonsingular local decay estimate one gets 
a pointwise bound for this integral that is at most quadratically growing in x and y and 
decaying like |t| _1,/2 in the even case and |t|~ 3//2 in the odd case. Since the estimates involve 
up to two derivatives of the quotient in the integrand, the bounds will be uniform in \i for 

To prove the decay for large r = \t\ in the truly singular case when \i > /i min > 0, we split 
k(x, y; t) into three parts. Let <?a(A) and <7a(A) be as before, and introduce the new "bump" 
functions gc(X) and gc(X) = 1 — <7g(A), both infinitely differentiable and nonnegative, with 
<?g(A) identically equal to zero outside of the interval (-JJI— 3G/2, ■ s /JI+3G/2) and identically 
equal to one inside of the interval (y/Ji — G/2, JJL + G/2). Set 

{x)-°{y)-< 



k{x,y\t) 



where 



lim 

54.0 



3A / 2 ^(x,0,A)^ o) (y,0,A)- 
A 2 -/i 

^°) {xAX )¥r\y,0,Xy 



■ 9A (X)e- 2lXt dX, 



X 2 



ik8 



9a{X)e 



-2iX 2 t 



dX 



(B.57) 



(B.58) 



(B.59) 



and 



v^/ 2 ^e. o >(x,0,A)^ o) (y,0,A)* 2ixH 
^ 2 _ gA{X)g G {A)e 



dX 



A/2 



+ 



^ o) (x > 0,A)^ o) (y,0 > A)* ^ 

9c{x)e 



;b.6o) 



dX 
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Note that in keeping the contributions near zero and near // distinct, we are assuming 
that 3 A/2 < ^//imin — G. The analysis of Iq and / proceeds exactly as in the proof of the 
nonsingular local decay estimate. The results are almost identical. For I one can integrate by 
parts as many times as one likes, and therefore one gets a pointwise estimate with arbitrary 
decay in time of order 0(|t|~ fe ) for k > 2, but at the cost of polynomial growth in x and y of 
degree k. For I , one gets a pointwise estimate that decays like 0(|t|~" 1//2 ) in the even case 
and 0(|t|" 3/2 ) in the odd case, at the cost of quadratic growth in x and y. 

The estimates of Iq and / are uniform for large /i. The pointwise bounds for I involve 
supremum bounds over the range of integration of the quantity 

A 2 ( Ak (Y^) ( A ) > (B.61) 

with /(A) given by ^'°\x, 0, \)¥°'°\y, 0, A)*# G (A) for A > JJL + G/2 and with /(A) given 
by ¥°'°\x,0, \)¥°' o) (y,0, \)*g A (\)g G (\) for A/2 < A < ^ - G/2. In particular, we will 
need these bounds for k = 2, in which case 

A2 U^) (A) = ( ^-T': 3 " ) fW+( :JS + 3 £ ) WHS®- (B-62) 



) 2 -Ai/ V A 2 (A 2 -^) 3 / VA(A 2 -Ai) 2 / A 2 -|i 

For the part of / involving A > it is easy to check that the three coefficients above 

are monotonic functions of A for A > y/Ji that decay for large A with \x fixed. Therefore each 
coefficient is bounded by its magnitude at the lower endpoint A = ^JJi J ^ G/2. With G held 
fixed, these bounds are then seen to be decaying functions of \i. For the part of / involving 
A G (A/2, y/Ji — G/2), it is easy to see that the coefficients blow up at both endpoints. 
Therefore, for A and G sufficiently small but independent of /i, the coefficients will be 
bounded by the maximum of their values at A = A/2 and A = yffi — G/2. Again, holding 
A and G fixed, one sees that the bounds are uniform for large /i. This direct argument 
shows that, at least for k = 2, the pointwise bound for / is uniform in /i. Establishing the 
uniformity of the pointwise estimate for I is easier; the denominator \ 2 — fi plays no essential 
role for A < 3 A/2 for \x sufficiently large. 

The new term that must be handled differently is 1^. 

Lemma B.12 For all k >2, the integral 1^ satisfies the pointwise estimate 

2 P k (x,y) 

^ k - 1 iip ' (R63) 

where Pk(x,y) is a polynomial in \x\ and \y\ of degree k with positive coefficients that are 
uniform in /i. 
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Proof: Consider first t > 0. Then, the quantity to estimate is 



I, = ljm/ ^ d ^' j^liv gG (x)e-"**d\ 



= 2ze- 2 ^lime m / ^' o) {x,0, \)¥°'°\y,0, \)*g G {\) / e" 2i(A2 -^ )s 

^° J^/Jl-ZG/2 Jt 

poo [ p^/Jl+3G/2 

= 2ie- 21 ^ / e 2 ^ / ^' o) {x,0,X)^ < d ,O \y,0,X)*g G {X)e- 2lX2s dX ds . 

Jt J^Ji-ZG/2 

(B.64) 

Now, with gc{X) vanishing to all orders at the integration endpoints, it is possible to bound 
the integral in square brackets by iterated integration by parts. The bound is 0(|s| _fc ) and 
grows in x and y like a polynomial Pk{x,y) of degree k. This bound is uniform in /i, since 
the only place /x appears is in the range of integration over which bounds are required, and 



from Lemma |B.2| we know that these bounds are uniform for all A. Therefore we have 

Pk(x, y) 
k-i It 



I M |<2/ P u (x,y)\8\- k d8= * 1 , (B.65) 



which establishes the lemma for t > 0. For t < 0, one gets an integral from — oo to t in the 
second step above, and ultimately obtains the same bound. □ 

Finally, we put the pieces together to complete the proof of Proposition p.2j . For the odd 



case, we want decay of order 0(|t|~ 3 / 2 ). For (y)~°~I to be in L 2 (M 2 ) with this decay 
rate, we need a > 5/2. With this bound on a, we can get (x)~°~ (y)~ a I being in L 2 (R 2 ) 
with decay bounded by 0(\t\~ 2 ) = o(|t|~ 3 / 2 ), but no better. Finally, for (x) _0 "(y) _0 "/ M to 
be in L 2 (R 2 ) with decay 0(|t|~ 2 ) we now see that we need to localize a bit more in space 
by taking a > 7/2. Combining these large time estimates with the finite time bound of 
Lemma [B.ll| establishes the proposition in the odd case. Similar arguments for the even 



case give an L 2 (IR 2 ) norm that decays like 0(|t| -1 / 2 ) for a > 7/2. This finishes the proof of 
Proposition |B.2j . 

Remark: Evidently, the singular decay estimates blow up when /x m i n approaches zero. This 
is an essential phenemenon in both the odd and even cases. This is best seen by considering 
the singular integral for the case fi = 0: 



^ O \x,o,x)^r\y^,xy _ 2lXH 







A 2 — ikS 



dX. (B.66) 



This integral blows up for all x, y, and t, as 6 tends to zero in the even case. In the odd case 
there is sufficient vanishing at A = for the limit of S I to exist for all x, y, and t, but the 
limit only decays in t like |t| -1 ^ 2 . Thus it is not possible for estimates of the form derived 
for > /i min > to hold uniformly in any neighborhood of \i = 0.D 
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